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Abstract 


Scheduling  can  be  fonnalized  as  a  constraint  satisfaction  problem  (CSP)‘  Within  this  framework 
activities  in  a  plan  are  interconnected  via  temporal  relation  constraints  a  la  Allen,  thereby 
defining  a  temporal  constraint  graph  (TCG).  Additionally  there  are  capacity  constraints 
restricting  the  use  of  each  resource  to  only  one  activity  at  a  time.  Together  these  constraints  form 
a  tempor all  capacity  constraint  graph  (T/CCG)  .  Preferences  such  as  meeting  due  dates, 
reducing  order  flowtime,  or  selecting  accurate  machines  are  modeled  as  utility  functions  over  the 
domain  of  p>ossible  start  times  and  durations  of  activities  and  over  the  sets  of  possible  resources 
activities  can  use.  These  preferences  interact  via  the  TCG  and  via  the  resource  edacity 
constraints.  Hence,  in  general,  they  cannot  be  simultaneously  optimized.  The  objeaive  of 
preference  propagation  techniques  is  to  transform  such  local  a  priori  preferences  so  as  to 
account  for  their  interactions. 

In  this  paper  we  describe  a  probabilistic  framework  in  which  start  time,  duration  and  resource 
preferences  are  propagated  across  T/CCGs  in  order  to  focus  attention  in  an  incremental 
scheduler.  The  propagation  is  first  performed  across  the  TCG,  thereby  producing  activity  (a 
posteriori)  start  time  and  duration  distributions.  These  distributions  allow  for  early  detection  of 
unsolvable  CSPs,  and  provide  measures  of  value  goodness  and  variable  looseness  for  aaivity 
start  times  and  durations.  In  a  second  phase,  these  distributions  are  combined  to  predia  t  z 
degree  of  contention  for  each  resource  and  the  reliance  of  each  activity  on  the  possession  of  that 
resource. 
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1.  Introduction 
1.1.  The  Issue 

We  are  concerned  with  the  issue  of  how  to  opportunistically  focus  an  incremental  scheduler’s 
attention  on  the  most  critical  decision  points  and  the  most  promising  decisions  in  order  to  reduce 
search  and  improve  the  quality  of  the  resulting  schedule.  More  specifically  we  are  concerned 
with  incremental  constraint  directed  scheduling  where  the  pros^iem  is  defined  as  a  set  of 
variables  and  a  set  of  constraints.  Both  variables  and  constraints  are  determined  by  the  initial 
scheduling  problem  and  the  earlier  decisions  made  by  the  scheduler.  The  interactions  of  the 
constraints  determine  the  structure  of  the  problem  space.  We  characterize  the  problem  space  with 
a  set  of  texture  measures  that  are  used  to  both  identify  critical  decision  points  and  selea  a 
decision  at  each  of  these  points.  The  process  of  analyzing  the  current  problem  and  generating 
new  decisions  (e.g.  scheduling  an  operation)  is  repeated,  thereby  resulting  in  the  incremental 
construction  of  a  schedule. 

Real-life  scheduling  problems  are  subject  to  a  variety  of  preferences  [Johnson  74,  Fox  83,  Ow 
84,  Smith  86]  such  as  meeting  due  dates,  reducing  the  number  of  machine  set-ups,  reducing 
inventory  costs,  using  accurate  and/or  fast  machines,  making  sure  that  some  jobs  are  performed 
within  a  single  work-shift,  etc.  Although  these  preferences  are  usually  set  independently  to  one 
another,  they  interact.  For  instance  selection  of  a  good  start  time  for  an  activity  (e.g.  to  meet  a 
due  date)  may  prevent  the  selection  of  an  accurate  machine  for  another  operation  or  may  prevent 
meeting  another  job’s  due  date.  For  this  reason,  selecting  operation  stan  times  or  allocating 
resources  based  solely  on  local  a  priori  preferences  is  likely  to  result  in  poor  schedules. 
Preference  propagation  is  meant  to  allow  for  the  constmction  of  measures  that  reflect  preference 
interactions.  These  measures  can  then  serve  to  guide  the  construction  of  a  good  overall  schedule 
rather  than  a  schedule  that  locally  optimizes  a  subset  of  preferences. 

We  perform  preference  propagation  within  a  probabilistic  fiamework.  We  associate  with  each 
variable’s  value  a  probability  that  reflects  the  likelihood  that  this  value  results  in  a  good  schedule 
overall  {value  goodness).  These  probabilities  are  refined  by  being  propagated  across  the  problem 
constraints.  Value  goodness  is  a  texture  measure  that  helps  selecting  assigrunents  for  variables. 
Identification  of  critical  variables  (i.e.  decision  points)  is  performed  using  another  texture 
measure,  called  variable  looseness.  A  critical  variable  or  group  of  variables  is  one  whose  good 
overall  values  are  likely  to  become  unavailable  if  one  were  to  start  assigning  values  to  other 
variables  first.  Notice  that  if  our  measures  of  value  goodness  were  perfect,  the  order  in  which 
variables  are  instantiated  would  not  matter.  However  such  perfect  measures  could  only  be 
obtained  by  first  solving  the  problem.  Because  in  practice  measures  of  value  goodness  contain 
some  uncertainty,  one  has  to  account  for  the  effects  of  assigning  a  value  to  a  variable  over  the 
availability  of  good  values  for  other  variables.  A  variable  instantiation  order  is  accordingly 
defined  starting  with  the  most  critical  (i.e.  least  loose)  variables.  In  this  paper  we  are  interested 
in  the  identification  of  critical  activities  (i.e.  operations).  An  activity  is  made  of  a  start  time 
variable,  possibly  a  duration  variable,  and  a  set  of  resource  variables.  We  identify  critical 
activities  as  the  ones  that  heavily  rely  on  the  possession  of  highly  contended  resources.  Indeed,  if 
such  critical  activities  are  not  scheduled  first,  it  is  very  likely,  that  by  the  time  the  scheduler  turns 
its  attention  to  these  activities,  the  resources  that  would  have  been  the  most  appropriate  for  these 
activities  will  no  longer  be  available. 
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We  discuss  preference  propagation  in  temporal/capacity  constraint  graphs  (T/CCG).  Temporal 
constraints  define  partial  orderings  among  the  activities  to  be  scheduled.  All  thirteen  of  Allen's 
[Allen  84]  temporal  relation  constraints  are  accounted  for.  Resource  capacity  constraints  restnet 
the  use  of  resources  to  only  one  activity  at  a  time.  Both  situations  with  fixed  and  variable 
duration  activities  are  discussed.  Our  formalism  allows  for  both  activity  start  time  and  duration 
preferences  as  well  as  for  resource  preferences.  It  also  accounts  for  prior  resource  reservations  if 
any.  It  is  shown  that  tht:  (a  posteriori)  start  time  and  Juration  distributions  resulting  from  the 
propagation  across  the  temporal  constraints  can  be  combined  to  identify  resources  that  are  highly 
contended  for  {resource  conteniion)  and  activities  that  heavily  rely  on  the  possession  of  these 
resources  {activity  resource  reliance)  in  function  of  time. 

We  also  argue  that  a  posteriori  start  time/duration  distributions  can  be  seen  locally  as  measures 
of  start  time/duration  goodness  and  globally  as  measures  of  start  time/duration  looseness.  Our 
notion  start  time/duration  looseness  generalizes  the  Operations  Research  notion  of  slack  [Baker 
74,  Johnson  74]. 


1.2.  Formalization  of  the  Scheduling  Problem 

The  faaory  scheduling  problem  is  often  described  as  a  two  step  problem:  a  process  planning 
step  and  a  resource  planning  step  [Fox  83].  Process  planning  deals  with  the  generation  and 
selection  of  plans  (i.e.  process  routings)  that  satisfy  the  order  specifications.  Resource  planning, 
sometimes  also  referred  to  as  scheduling,  deals  with  the  allocation  of  resources  (e.g.  machines) 
to  activities  and  the  assignment  of  start  and  end  times  to  activities.  In  general  both  steps  can  be 
interleaved. 

In  this  paper  we  will  be  concerned  exclusively  with  the  scheduling  part  of  the  problem:  we  wlU 
assume  that  we  are  given  a  set  of  plans  to  schedule.  Here  a  plan  is  simply  defined  as  a  partial 
ordering  of  activities.  Each  activity  may  require  one  or  more  resources,  for  each  of  which  there 
can  be  several  alternatives. 

We  formalize  the  scheduling  problem  as  a  constraint  satisfaction  problem  (CSP).  The  variables 
of  the  problem  are  the  activity  start  times,  the  resources  allocated  to  each  activity,  when  there  is  a 
choice,  and  possibly  the  duration  of  each  activity.  An  activity’s  end  time  is  defined  as  the  sum  of 
the  activity’s  start  time  and  duration.  We  differentiate  between  two  types  of  constraints:  required 
constraints  and  preferential  constraints  [Fox  83].  Required  constraints  determine  the 
admissibility  of  a  solution  to  the  CSP  (schedule)  while  preferential  constraints  allow  for 
differentiating  among  admissible  solutions.  The  degree  of  satisfaction  of  a  preferential 
constraint  is  defined  by  a  utility  function  that  maps  the  possible  values  of  a  variable  onto  utilities 
ranging  between  0  and  1.  A  utility  ot  0  indicates  a  non-admissible  value.  A  value  with  utility  1  is 
an  optimal  value. 

We  will  be  dealing  explicitly  with  two  types  of  required  constraints:  temporal  relation 
constraints  and  resource  capacity  constraints.  Temporal  relation  constraints  are  used  to  describe 
partial  orderings  among  activities  as  provided  by  the  process  planning  step.  We  will  be  using 
Allen’s  temporal  relation  constraints  [Allen  84]  to  describe  these  constraints  (Figure  1-1).  We 
will  refer  to  the  graph  defined  by  these  constraints,  for  a  given  CSP,  as  the  CSP’s  temporal 
constraint  graph  (TCG).  Capacity  constraints  restrict  the  number  of  reservations  of  a  resource 
over  any  time  interval  to  the  capacity  of  that  resource.  In  this  paper,  for  the  sake  of  simpliciw. 
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Pictorial  Representation 


X  BEFORE  Y  XXX  YYY 

Y  AFTER  X 


X  EQUALS  Y 
Y  EQUALS  X 


XXX 

YYY 


X  MEETS  Y  XXXYYY 

Y  MET -BY  X 


X  OVERLAPS  Y  XXX 

Y  OVERLAPPED -BY  X  YYY 

X  DURING  Y  XXX 

Y  CONTAINS  X  YYYYY 

X  STARTS  Y  XXX 

Y  STARTED-BY  X  YYYYY 


X  FINISHES  Y  XXX 

Y  FINISHED-BY  X  YYYYY 


Figure  1-1:  Allen’s  13  temporal  relation  constraints 

we  will  always  be  assuming  resources  with  unary  capacity.  Together  these  required  constraints 
form  a  temporal/c2q)acity  constraint  graph  (T/CCG).  A  schedule  that  does  not  satisfy  the  required 
constraints  of  the  CSP  is  not  admissible. 

We  will  allow  for  preferential  constraints  on  activity  start  times  and  durations  as  well  as  on  the 
resources  to  be  used  by  each  activity.  Preferential  constraints  are  described  with  utility  functions. 
For  a  given  preferential  constraint,  a  variable’s  value  is  admissible  only  if  its  utility  is  strictly 
positive.  High  preference  for  an  admissible  value  is  indicated  by  a  high  utility.  In  practice  the 
domain  of  admissible  values  resulting  from  these  preferential  constraints,  i.e.  the  domain  with 
strictly  positive  utilities,  is  always  bounded.  For  instance  the  domain  of  admissible  start  times  of 
an  activity  is  constrained  at  one  end  by  the  order  release  date  and  at  the  other  end  by  the  order 
due  date  according  to  the  durations  of  the  activities  that  precede/follow  the  activity  within  the 
plan. 

Notations 

We  have  to  schedule  a  set  of  activities  {A,,  Aj,...,  A,}.  Let  denote  the  time  interval  over 
which  A^  spans,  st^,  et*  and  du^  respectively  denote  I^’s  start  time,  end  time,  and  duration. 
Activities  are  connected  by  a  set  of  temporal  relation  constraints,  thereby  forming  a  TCG.  We 
view  TCGs  as  undirected  graphs.  An  arc  in  a  TCG  indicates  the  presence  of  a  temporal  relation 

between  two  intervals  (e.g.  L  BEFORE  L  or  equivalently  Ij  AFTER  I,).  Let  C,,  Cj . denote 

the  temporal  relation  constraints  in  the  TCG.  The  TCG,  which  has  been  produced  during  the 
process  planning  phase,  is  assumed  consistent. 
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Additionally  there  are  capacity  constraints  limiting  the  use  of  each  resource  to  only  one 
activity  at  a  time.  By  adding  these  capacity  constraints  to  the  TCG,  one  obtains  the  CSP’s 
T/CCG.  In  this  paper  we  will  not  need  to  formalize  the  description  of  T/CCGs  any  further  as  we 
will  be  mainly  dealing  with  their  temporal  abstractions,  i.e.  the  TCGs  obtained  by  omitting  the 
capacity  constraints  in  the  T/CCGs. 

Each  activity  has  a  preferential  start  time  constraint  with  associated  utility  function  denoted 
u^^.  Activity  A^^’s  duration  is  either  fixed  or  constrained  by  a  preferential  constraint  with  utility 

function  u^  .  The  ranges  of  admissible  start  times  and  durations  are  assumed  to  be  bounded, 

k 

which  is  always  the  case  in  practice. 

Each  activity  A^  may  require  one  or  more  resources  R^,,  R^.  For  each  resource  R^^ 

required  by  activity  A^,  there  is  a  set  of  possible  resources  R^j,  available  on  the 

factory  floor.  This  set  is  assumed  to  be  finite,  which  is  also  always  the  case  in  practice.  These 
different  resources  are  usually  not  equally  preferred.  A  resource  utility  function,  u^^  associates  a 
utility  (preference)  ngJR^^p  with  each  possible  resource  R^^y.  For  instance  a  milling  operation 
may  require  a  milling  machine  and  a  human  operator.  There  may  be  two  milling  machines 
available  on  the  factory  floor.  For  this  specific  milling  operation,  nulling-machine  j  may  have  a 
preference  of  1.0  and  milling  machine2  a  preference  of  0.4  (e.g.  due  to  a  difference  in  the 
accuracy  of  the  machines.).  There  may  dso  be  several  human  operators  available  with  different 
utilities. 

The  global  utility  of  a  schedule  is  obtained  by  summing  all  the  preferential  constraints  ’  utilities 
(for  the  given  schedule). 

U.  An  Example 

We  now  introduce  a  simple  scheduling  problem  that  we  will  use  throughout  this  p^r. 

The  problem  involves  scheduling  two  orders:  orderj  and  order^  (Figure  1  -2); 

•  order,  comprises  five  activities:  Aj,  A2,...rA.^, 

•  order2  comprises  three  activities:  Ag,  Aj,  and  Ag. 

All  activities  have  the  same  duration,  namely  30  time  units.  Ci,C2>—  C7  are  the  temporal 
relation  constraints  imposed  by  the  process  planning  step.  For  instance,  indicates  that  A^  has 
to  precede  A2.  The  domain  comprises  three  physical  resources:  Rj,  R2,  and  R3. 


•  Aj  requires  a  resource  Rjj  which  can  be  either  Rj  or  R2  (with  equal  preference),  i.e. 

anduj^^^(R3)=0. 


•  A2  requires  a  resource 

^21 

which 

has 

to 

be 

Rj,  i.e. 

Ur^^(Ri)=1, 

and 

•  A3  requires  a  resource 

R31 

which 

has 

to 

be 

R2,  i.e. 

and 

“R3,(^i)='^R3,(^3)=^- 

•  A4  requires  a  resource 

R41 

which 

has 

to 

be 

R2,  i.e. 

and 

Ur^^(Ri)=Ur^_(R3)=0. 
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•  Aj  requires  a  resource  R^j  which  has  to  be  Rj,  i.c.  Uj^  (R|)=l,  and 

•  requires  a  resource  R^j  which  has  to  be  R3,  i.e.  Uj^  (R3)=l,  and 

•  A-j  requires  a  resource  R-^j  which  can  be  either  R2  or  R3  (with  equal  preference),  i.e. 

“R7,(R2)=“R7,(R3)=1-  “R7,(Ri)=<^- 

•  Ag  requires  a  resource  Rgj  which  has  to  be  R3,  Le.  Uj^  (R3)=l,  and 

“Rg,(^l)=“Rg,(^2)=^- 


C,:  A,  BEFORE  Az 
G:  A 2 BEFORE  A, 
Cji  A  3  BEFORE  A  5 
C4:  A,  BEFORE  A* 
C5:  A  4  BEFORE  Aj 
Q:  A «  BEFORE  A7 
C,:  A  7  BEFORE  A« 


Figure  1-2:  TCG  for  a  two  order  scheduling  problem 

Activities  in  order ^  (Le.  A^,...,  A5)  are  assumed  to  have  the  same  start  ime  utility  function. 
TTie  function  requires  that  these  activities  start  between  time  0  and  time  140,  with  an  optimal 
start  time  at  120  (Figure  1-3).  Time  0  can  be  interpreted  as  the  order  release  date.  Time  140  +  30 
=  170  (latest-start-time  +  duration  =  latest-finish-time)  could  represent  the  time  after  which  the 
client  would  refuse  the  order. 

,MI  activities  in  order2  will  be  assumed  to  have  a  unifomi  start  time  utility  between  0  and  120 
igure  1-3). 


general  preferential  constraints  are  set  independently  to  one  another  and  may  therefore  be 
incompatible.  For  instance,  it  will  obviously  be  impossible  to  simultaneously  schedule  both  A-, 
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I 


0 


Figure  1-3:  Start  time  utility  functions 

and  A3  at  their  optimal  start  time,  namely  120.  Therefore  instead  of  a  priori  preferences,  one 
needs  (a  posteriori)  preferences  indicating  values  that  are  likely  to  result  in  a  good  schedule 
overall.  Such  preferences  can  only  be  obtained  by  accounting  for  constraint  interactions.  This  is 
the  objective  of  the  propagation  technique  presented  in  this  paper. 

Looking  more  closely  at  our  example,  one  notices  that  there  are  four  activities  requiring 
resource  R2:  Aj,  A3,  A4,  and  A^.  Out  of  these  four  activities,  up  to  three  can  occur  in  parallel, 
namely  A3,  A4,  and  A^.  These  activities  will  therefore  compete  for  the  possession  of  R2-  There  is 
no  such  competition  for  R^  and  R3  as  the  activities  that  require  these  resources  are  fully  ordered 
temporally  (e.g.  Aj  has  to  be  carried  out  before  A2).  The  scheduling  of  R2  is  therefore  more 
critical  than  that  of  Rj  and  R3.  An  incremental  scheduler  should  first  focus  its  attention  on  the 
competition  between  A3,  A4,  and  A7  for  R2.  Moreover,  since  A7  has  two  resource  alternatives 
(R2  and  R3)  while  A3  and  A4  have  only  one,  and  since  A3  has  less  slack  than  A4,  we  would  like 
our  incremental  scheduler  to  first  schedule  A3  with  R2^. 

A  more  detailed  analysis  confirms  that  first  scheduling  A3  rather  than  A4  is  the  right  decision. 
It  also  reveals  the  influence  of  the  preferential  constraints  (in  this  case  the  start  time  preferences) 
in  determining  activity  criticality.  We  consider  two  scenarios:  one  where  the  first  activity  to  be 
scheduled  is  A4  (scenariol)  and  one  where  it  is  A3  (scenarioZ).  Given  that  A5  cannot  start  later 
than  140  and  that  A3  and  A4  have  a  duration  of  30,  A^,  A2,  A3,  and  A4  cannot  start  later  than 
1 10.  Hence,  according  to  their  start  time  utility  functions,  these  activities  will  prefer  to  start  as 
late  as  possible  (Figure  1-3).  In  scenariol,  A4  is  the  fiist  activity  to  be  scheduled.  It  is  scheduled 
as  late  as  possible  (while  still  leaving  some  room  for  A5  to  have  a  good  schedule),  say  at  time  90. 
Since  both  A3  and  A4  require  R2,  A3  has  to  be  schedule  before  A4.  The  resulting  schedule  is  the 
sequence  A^,  A2,  A3,  A4,  and  A5,  as  displayed  in  Figure  1-4.  Alternatively  (scenarioZ)  suppose 
that  we  decide  to  first  schedule  A3.  This  time  A3  is  scheduled  to  start  at  time  90,  and  A4  has  to 
occur  before  A3.  However  and  A4  can  occur  in  parallel  before  A3  (Figure  1-5).  Hence, 
globally,  Aj  to  A4  start  later  in  scenarioZ  than  in  scenariol.  In  other  words,  scenarioZ  results  in 
a  higher  global  utility  than  scenariol. 

'Notice  that  we  arc  assuming  an  incremental  scheduler  whose  reservations  are  nonpreempuble.  The  order  in 
which  activities  are  allocated  resources  would  not  maner  if  allocations  were  preemptible.  Most  predictive 
schedulers  do  not  allow  for  such  preemptions  as  they  tend  to  produce  infinite  loops  if  one  docs  not  take  special 
precaution. 
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Our  approacli  to  preference  propagation  formalizes  the  above  considerations. 
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Figure  1-4:  Gantt  chart  for  a  schedule  of  orderl  obtained  with  scenariol. 

The  global  start  time  utility  has  been  obtained  by  adding  the  start  time 
utility  of  the  five  activities. 
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Figure  1-5:  Gantt  chart  for  a  schedule  of  orderl  obtained  with  scenario2. 

The  global  start  time  utility  has  been  obtained  by  adding  the  start  time 
utility  of  the  five  activities. 


1.4.  Organization  of  the  Paper 

In  the  next  section  we  give  an  overview  of  related  work  in  constraint  satisfaction  and 
scheduling.  Section  3  introduces  the  assumptions  that  arc  the  basis  to  our  probabilistic  approach 
to  preference  propagation  and  gives  an  overview  of  the  approach.  Section  4  describes  the 
propagation  of  start  time  and  duration  probability  distributions  in  TCGs.  As  already  mentioned 
earlier,  when  one  is  simply  concerned  with  knowing  whether  two  time  intervals  are  temporally 
related  or  not,  a  TCG  can  be  considered  as  an  undirected  graph.  When  we  will  be  talking  about 
cycles  in  the  TCG.  we  will  always  be  referring  to  the  undirected  interpretation  of  the  graph. 
Subsection  4.2  deals  with  acyclic  TCGs  with  fixed-duration  activities.  Subsection  4.3  relaxes 
the  fixed-duration  hypotliesis.  Subsection  4.4  relaxes  the  acyclicity  assumption.  Subsection  4.5 
discus.scs  propagation  in  general  TCGs  where  tliere  are  explicit  disjunctions  of  temporal  relation 
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constraints.  Finally  subsection  4.6  analyzes  the  results  of  the  section  with  respect  to  a  set  of 
requirements  and  desiderata  identified  in  subsection  3.2.  Section  5  explains  how  the  results  of 
the  previous  section  can  be  combined  to  obtain  measures  of  resource  contention  and  activity 
resource  reliance.  In  section  6,  we  discuss  the  time  complexity  and  the  expressiveness  of  our 
framework  as  well  as  possible  improvements.  Section  7  summarizes  the  main  ideas  of  the  paper. 

For  the  sake  of  concision,  subsections  4.2  and  4.3  contain  only  the  treatment  of  two  temporal 
relation  constraints  (BEFORE,  and  MEETS).  Formulas  for  the  complete  set  of  Allen’s  temporal 
relation  constraints  is  presented  in  appendix  1.  For  the  same  reason,  section  5  only  sketches  the 
computation  of  resource  demand  densities.  The  reader  will  find  a  complete  treatment  of  these 
densities  in  appendix  2. 

2.  Related  Work 

We  already  mentioned  that  this  paper  does  not  deal  with  process  planning.  Hence  it  is 
assumed  that  the  TCGs  to  be  scheduled  are  consistent.  [Vilain  86]  has  proved  that  consistency 
checking  in  a  general  TCG  is  NP-hard.  Several  algorithms  have  been  proposed  in  the  literature  to 
perform  partial  or  total  consistency  checking  in  a  TCG.  Allen’s  algorithm  [Allen  83]  achieves 
3 -consistency^  in  a  general  TCG  in  polynomial  time  and  space.  A  complete  consistency 
checking  algorithm  using  a  variation  of  data  dependency  backtracking  is  presented  in 
[Valdes-Perez  87].  Although  the  algorithm  is  designed  for  quick  pruning,  its  asymptotic 
complexity  remains  exponential.  [Vilain  86]  and  [Tsang  87]  point  out  that  consistency  checking 
can  actually  be  perfoimed  in  polynomial  time  provided  that  the  TCG  does  not  contain  cenain 
types  of  disjunctive  relations  such  as  "Intervalj  has  to  be  either  BEFORE  or  AFTER  Interval,". 

When  additional  constraints  such  as  capacity  constraints,  preferential  start  time  constraints, 
preferential  duration  constraints,  preferential  resource  constraints,  or  resource  reservations  are 
added  to  a  consistent  TCG,  the  resulting  CSP  may  stop  being  consistent.  These  are  the  types  of 
inconsistencies  that  we  will  be  referring  to  later  in  this  paper.  Consider  the  simple  example 
depicted  in  Figure  2-1.  There  are  two  activities:  Aj  and  A^.  Aj  is  BEFORE  A^.  Aj  has  a 
preferential  start  time  constraint  specifying  that  it  has  to  start  between  10  and  15.  Aj’s 
preferential  start  time  constraint  specifies  that  Aj  has  to  start  between  0  and  5.  A,’s  duration  is 
10.  Aj’s  duration  does  not  matter.  The  resulting  CSP  is  obviously  inconsistent  (unsatisfiable)  as 
Aj’s  earliest  end  time  (10+10=20)  is  after  Aj’s  latest  start  time  (5). 

Propagation  of  activity  start  and  end  time  windows  (earliest/latest  start  and  end  times)  dates 
back  to  the  CPM  algorithm  [Johnson  74].  The  PERT  method  generalizes  CPM  by  allowing  for 
uncertainty  in  activity  durations.  [Vere  83]  adapted  the  CPM  propagation  techniques  to  a 
planning  system  called  DEVISER.  In  DEVISER  start  time  windows  are  described  as  triples  of 
the  form  {earliest-start-time,  ideal-start-time,  latest-start-time).  The  ideal  start  time  information 
is  optional.  A  start  time  triple  can  be  seen  as  a  triangle-shaped^  start  time  utility  function  (Figure 
2-2).  When  the  ideal  start  time  is  omitted,  the  window  can  be  interpreted  as  a  rectangle-shaped 


^According  to  [Freuder  82],  a  constraint  graph  is  k-consistent  if  for  any  set  of  (k-1)  variables,  any  consistent 
assignment  of  values  to  these  (k-1)  variables,  and  any  k-th  variable,  there  always  exists  a  value  for  the  k-th  variable 
such  that  the  k  values  taken  together  (i.e.  for  the  (k-l)+l  variables)  are  consistent. 

^This  is  not  the  only  possible  interpretation. 


C:  Ai  BEFORE  A2 
du  , :  A,  's  duration 

Figure  2-1:  An  inconsistent  CSP  with  a  consistent  TCG 

start  time  utility  function.  In  DEVISER  stan  time  windows  are  dynamically  compressed  to 
account  for  new  temporal  relations  and  activities  introduced  during  the  planning  process.  The 
compression  does  not  account  for  ideal  start  times  which  remain  fixed.  DEVISER  is  able  to 
handle  both  BEFORE/AFTER  and  MEETS/MET-BY  relations.  Because  preferences  are  not 
propagated,  the  system  performs  poorly  when  it  has  to  select  a  start  time  within  a  compressed 
window:  the  selection  is  based  on  purely  local  information. 


case  widi  an  ideal  ca.se  withoui  ideal 

start  time  start  time 


time  time 


est:  earliest  start  time 
1st:  latest  start  time 
ist:  ideal  start  time 


Figure  2-2:  A  utility  interpretation  of  DEVISER’S  st;ut  time  windows 

A  variation  of  Vere’s  algorithm  is  presented  in  [Bell  84]  that  accounts  for  viu-iable  duration 
activities.  (Smitii  83]  extends  Vere’s  approach  by  :iiso  accounting  for  DURING/CONTAl.MS 
temporal  relation  constraints  (see  also  fI.^Pape  87|),  Smith’s  temporal  module  also  ;iJlows  for  tlic 
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propagation  of  resource  reservations  over  T/CCGs.  [Rit  86]  describes  a  Waltz  algorithm  to 
propagate  generalized  temporal  windows  over  TCGs.  A  generalized  window  is  the  composition 
of  a. start  time,  end  time,  and  duration  window  (Figure  2-3).  The  method  accounts  for  all  13  of 
Allen’s  temporal  relations  as  well  as  for  disjunctions  among  these  relations.  In  the  general  case, 
Rit’s  algorithm  is  only  guaranteed  to  achieve  arc -consistency  [Mackworth  77].  Total  consistency 
can  however  be  guaranteed  in  TCGs  that  do  not  contain  disjunctions  of  temporal  relation 
constraints. 


St:  start  time 
et:  end  time 
est:  earliest  start  time 
1st:  latest  start  time 
eet:  earliest  end  time 
let:  latest  end  time 
md:  minimum  duration 
Md:  maximum  duration 


0  est  1st  St 

Figure  2-3:  Rit’s  generalized  window 

None  of  the  propagation  techniques  that  we  just  described  handles  preferences.  In  praaice, 
however,  different  times  within  a  window  are  not  equally  preferred.  For  instance,  in  the  factory 
scheduling  domain,  due  dates  and  associated  late  delivery  penalties  induce  preferences  on 
activity  end  times  (and  hence  start  times).  Inventory  costs  are  another  source  of  start  time 
preferences.  In  general  preferences  cannot  be  accounted  for  independently.  Selection  of  the 
optimal  Stan  time  for  one  activity  may  prevent  selection  of  the  optimal  start  time  for  another 
either  because  of  temporal  relation  constraints  between  the  two  activities  or  because  of  capacity 
constraints,  or  a  combination  of  the  two.  This  is  why  it  is  crucial  not  only  to  propagate  time 
windows  but  also  to  propagate  preferences  over  these  windows.  Window  propagation  simply 
guarantees  admissibility  of  the  values  within  a  compressed  window.  Preference  propagation 
strives  not  only  for  admissibility  but  also  for  optimality  by  locally  reflecting  preference 
interactions. 

Our  purpose  is  to  develop  preference  propagation  techniques  to  guide  an  incremental 
scheduler.  Both  empirical  and  analytical  studies  reported  in  [Haralick  80,  Freuder  82,  Purdom 
83,  Nadel  86a,  Nadel  86b,  Nadel  86c,  Stone  86]  indicate  that,  in  general,  the  amount  of  search 
required  to  find  a  solution  to  a  CSP  can  be  significantly  reduced  by  using  the  following  two 
look-ahead  schemes  [Dechter  88]; 

1.  Variable  Ordering:  Focus  on  the  most  constrained  variables  first. 

2.  Value  Ordering:  Try  the  least  constraining  values  first. 

Tightly  constrained  variables  and  constraining  values  are  determined  by  the  interactions  of  the 
problem  constraints.  In  CSPs  where  variables  have  finite  sets  of  possible  values  and  all  values 
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are  equally  preferred,  the  number  of  possible  values  left  after  constraint  propagation  (i.e. 
consistency  checking)  can  be  used  to  determine  variable  tightnessAooseness.  In  problems  where 
values  are  not  equally  preferred,  like  in  the  scheduling  domain,  this  is  not  sufficient.  One  has 
also  to  account  for  value  goodness,  i.e.  the  utility  of  a  value  and  the  impact  of  selecting  that 
value  on  the  availability  of  good  values  for  the  other  variables.  Our  notion  of  value  goodness 
and  our  generalization  of  the  notion  of  variable  looseness  are  intended  to  allow  for  the 
generalization  of  these  two  look-ahead  strategies  to  CSPs  where  variables  can  have  infinite 
bounded  sets  of  possible  values  with  non-uniform  preferences. 

[Muscettola  87]  presents  a  probabilistic  framework  to  compute  resource  contention  in 
T/CCGs  with  only  BEFORE/APTER  temporal  relation  constraints  based  on  assumptions  on  the 
order  in  which  the  activities  are  scheduled.  This  paper  extends  Muscettola ’s  approach  for 
computing  resource  contention  by  removing  the  need  for  assumptions  on  the  order  in  which 
activities  are  scheduled,  by  dealing  with  ail  of  Allen’s  constraints,  by  allowing  for  duration  and 
resource  preferences  and  by  accounting  explicitly  for  earlier  resource  reservations. 

3.  A  Probabilistic  Framework  for  Preference  Propagation 

3.1.  An  Overview  of  the  Approach 

Our  purpose  is  to  develop  preference  propagation  techniques  to  guide  an  incremental 
scheduler.  An  incremental  scheduler  works  by  iterating  through  a  two-phase  process.  In  the  first 
phase  it  analyzes  the  structure  of  the  CSP  resulting  from  the  initial  scheduling  problem  and  the 
decisions  that  have  already  been  made.  In  the  second  phase,  based  on  this  analysis,  new 
decisions  are  generated  resulting  in  the  expansion  of  the  current  schedule  (e.g.  new  activities  are 
scheduled).  If  the  scheduler  reaches  a  deadend,  it  backtracks.  The  process  goes  on  until  a 
satisfactory  schedule  is  produced. 

The  first  phase  analysis  is  performed  using  preference  propagation  to  dynamically  identify 
critical  decision  points.  Such  decision  points  are  determined  by  the  interactions  of  the  problem 
constraints.  In  the  case  of  the  scheduling  problem,  there  are  two  main  types  of  interaaions: 
operation  precedence  interactions  and  resource  requirement  interactions  [Smith  85]. 

1,  The  operation  precedence  interactions  are  the  ones  induced  by  the  TCG.  They  are 
sometimes  also  referred  to  as  intra-order  interactions'^. 

2.  Resource  requirement  interactions  are  induced  by  the  capacity  constraints.  They 
arise  from  the  contention  of  several  activities  for  the  same  resource.  They  are 
sometimes  referred  to  as  inter-order  interactions. 

Intra-order  and  inter-order  interactions  have  respectively  motivated  so-called  order-based  and 
resource-based  scheduling  techniques.  In  the  past  few  years  it  has  become  clear  that  efficient 
scheduling  requires  the  ability  to  combine  these  two  perspectives  [Smith  85,  Smith  86]  so  as  two 
account  for  both  types  of  interactions.  Unfortimately  both  types  of  interactions  are  not  totally 
independent.  Intra-order  interactions  affect  the  time  intervals  over  which  activities  will  contend 
for  resources,  thereby  influencing  inter-order  interactions.  Resource  contention  in  turn  restricts 


‘‘Although  activities  within  a  same  order  can  also  interact  by  competing  for  the  same  resources. 
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the  times  over  which  activities  can  occur,  thereby  influencing  intra-order  interactions. 

In  order  to  deal  with  the  uncertainty  in  the  interactions  between  uninstantiated  variables,  we 
have  adopted  a  probabilistic  model.  For  each  uninstantiated  variable,  a  probability  density  is 
computed  for  the  variable’s  possible  values  that  indicates  the  likelihood  of  each  value  to  result  in 
a  good  schedule  overall,  given  the  decisions  already  made  by  the  incremental  scheduler.  This 
probability  is  a  dynamic  measure  of  value  goodness.  When  the  corresponding  probability  density 
is  normalized,  we  also  interpret  this  probability  as  the  dynamic  probability  that  the  scheduler 
assigns  that  value  to  the  vari^le. 

In  its  simplest  form  our  approach  involves  the  following  steps: 

1.  Based  on  the  current  partial  schedule  as  well  as  start  time,  duration  and  resource 
preferences,  a  priori  probability  distributions  are  produced  for  the  start  time, 
duration  and  resources  of  each  unscheduled  activity, 

2.  These  a  priori  probability  distributions  are  propagated  over  the  TCG,  resulting  in  a 
posteriori  start  time  and  duration  probability  distributions, 

3.  The  a  posteriori  distributions  obtained  in  the  previous  step  are  combined  to 
compute  activity  individual  demand  densities.  An  activity  A^’s  individual  demand 
density  at  time  t  for  a  resource  say  D^J(t),  is  defined  as  the  probability  that 

is  active  at  time  t  and  uses  R^j  to  fulfill  its  resource  requirement 

4.  Activity  individual  demand  densities  are  combined  to  measure  resource  aggregate 
demand  densities.  The  aggregate  demand  density  for  a  resource  at  time  t  is  the 
probabilistic  demand  for  that  resource  at  time  t. 

Both  iterative  and  hierarchical  variations  of  this  basic  propagation  algorithm  will  be  discussed. 

From  an  Operations  Research  point  of  view,  activity  a  posteriori  start  time  and  durations 
provide  a  measure  of  intra-order  interactions  and  a  generalization  of  the  notion  of  slack  [Johnson 
74].  Resource  aggregate  demand  densities  reflect  the  level  of  resource  contention  defined  by  the 
CSP’s  inter-order  interactions.  Resource  aggregate  demand  densities  can  be  identified  with  the 
Operations  Research  concept  of  bottleneck  analysis  [Smith  85,  Smith  86,  Muscettola  87]. 

From  a  constraint  satisfaction  perspective,  a  posteriori  start  time/duration  distributions  should 
be  regarded  locally  as  measures  of  value  goodness,  and  globally  as  measures  of  variable 
looseness  for  activity  start  times  and  durations.  Aggregate  demand  densities  can  be  interpreted  as 
measures  of  constraint  contention  and  individual  demand  densities  as  measures  of  activity 
resource  reliance. 


3.2.  Building  A  Priori  Probability  Distributions  Based  on  Local  A  Priori  Preferences 
Our  approach  uses  Bayesian  probabilities  to  estimate  value  goodness,  i.e.  the  likelihood  that  a 
given  value  will  result  in  a  good  schedule  overall.  It  consists  in  the  construction  of  a  priori 
probability  distributions  for  each  uninstantiated  variable  based  on  local  preferences.  These 
probabilities  are  then  refined  so  as  to  account  for  constraint  interactions,  thereby  resulting  into  a 
posteriori  probability  distributions. 

Obviously  the  main  concern  in  such  an  approach  is  to  obtain  good  estimates  of  value  goodness 
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(desideratuml)  as  these  estimates  are  essential  to  both  the  identification  of  critical  decision 
points  in  the  search  space  and  the  selection  of  a  decision  at  these  points.  There  are  however 
some  more  specific  requirements  and  desiderata  to  which  one  should  give  special  attention. 
Indeed,  besides  its  need  for  good  focus  of  attention  mechanisms,  efficient  search  also  requires 
the  ability  to  quickly  pmne  deadend  paths  in  the  search  tree. 

In  our  probabilistic  framework,  unsatisfi^iiity  is  detected  when  a  posteriori  probability 
densities  are  uniformly  zero.  This  indicates  that  the  interactions  of  the  problem  constraints  have 
reduced  the  set  of  admissible  values  for  a  variable  to  the  empty  set.  Detecting  unsatisfiability  in 
this  fashion  requires  that: 

1.  Requirement!:  Every  value  that  is  a  priori  admissible^is  given  a  strictly  positive  a 
priori  probability,  though  possibly  very  small. 

2.  Requirement2:  The  propagation  step,  which  combines  a  priori  probabilities  to 
account  for  constraint  interactions,  produces  a  posteriori  probability  densities  that 
are  zero  only  for  values  forbidden  by  the  constraint  interactions. 

Requirement!  restricts  the  construction  of  a  priori  probability  distributions.  Requirement2  is  a 
restriction  on  the  propagation  method  itself. 

Desideratuml:  Additionally,  in  order  to  detect  and  prune  inconsistent  states  as  soon  as 
possible,  one  would  like  a  posteriori  probabilities  to  be  zero  for  all  non-admissible  values 
(complete  consistency  checking).  Unfortunately  interactions  between  the  resource  requirements 
of  unscheduled  activities  seem  computationally  very  expensive  to  totally  account  for. 
Consequently  we  will  have  to  settle  for  partial  consistency  checking. 

These  general  requirements  and  desiderata  having  been  identified,  we  turn  our  attention  to  the 
constmction  of  a  priori  probability  distributions.  We  start  with  some  general  observations. 

In  the  presence  of  a  unique  variable  with  a  single  (unary)  preferential  constraint,  one  can  just 
select  one  of  the  optimal  values  defined  by  the  utility  function.  The  probability  density  for  the 
variable’  s  value  consists  of  a  set  of  peak  distributions  (Dirac  distributions),  each  centered 
around  one  of  the  optimal  values  (Figure  3- lb). 

On  the  other  hand,  in  the  presence  of  several  variables  and  constraints,  it  is  not  always  possible 
anymore  to  simultaneously  select  an  optimal  value  for  each  variable.  For  instance,  it  is  not 
always  possible  to  schedule  an  activity  at  its  optimal  start  time  and  with  a  set  of  optimal 
resources.  Very  often  one  has  settle  for  suboptimal  start  times  and/or  resources  in  order  to  find  a 
feasible  schedule  (i.e.  satisfy  all  the  CSP’s  constraints).  If,  for  a  given  variable,  the  interactions 
defined  by  the  constraints  are  weak,  it  is  usually  possible  to  select  a  value  that  is  stiU  very  close 
to  the  optimum  (or  one  of  the  optimums).  As  interactions  become  stronger,  it  becomes  more 
difficult  to  select  values  close  to  the  optimums:  the  probability  density  widens  (Figure  3-lc  and  d 
).  In  situations  of  extremely  strong  interactions,  one  is  just  happy  to  find  a  solution  within  the 
domain  of  admissibility  (non-zero  utility  value).  Hence  the  probability  distribution  tends  towards 
a  uniform  distribution  over  the  range  of  admissible  values  (Figxne  3-le). 


^Of  course,  a  more  sophisticated  analysis  will  result  in  the  rejection  of  a  larger  number  of  possible  values  (see 
desideratum2).  Therefore  what  is  really  important  for  requirement!  is  that  no  value  received  a  zero  a  priori 
probability  while  it  could  have  resulted  in  an  admissible  schedule. 
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Figure  3-1:  A  priori  probability  density  P(x)  for  a  variable  x  with  utility  function  u(x). 

This  is  an  example  with  a  single  optimal  value. 

The  a  priori  probability  densities  that  we  are  currently  using  are  essentially  obtained  by 
normalizing  utility  functions.  Intuitively  this  can  be  interpreted  as  a  sort  of  average  difficulty 
assumption  (see  Figure  3-1).  The  resulting  a  priori  probabilities  obviously  satisfy 
requirement!.  Moreover,  because  we  assume  that  domains  of  admissibility  defined  by  utility 
functions  are  bounded,  normalization  is  always  possible. 

Prior  to  normalizing  a  utility  function,  its  domain  is  pruned  to  account  for  earlier  resource 
reservations.  This  improves  the  quality  of  the  probability  distributions  with  respect  to 
desideratum2.  In  panicular  we  remove  from  the  start  time  probability  distributions  the  start 
times  that  are  not  allowed  by  the  current  resource  reservations.  A  start  time  t  is  not  allowed  for 
an  activity  A^.,  if  there  is  at  least  one  resource  required  by  A^  such  that  none  of  the  resources 
R. is  totally  available  between  t  and  r+5t  ,  where  5u  is  A.’s  smallest  admissible  duration. 

We  are  currently  investigating  alternative  methods  for  producing  a  priori  probability 
distributions.  In  particular  we  arc  investigating  both  iterative  and  hierarchical  approaches  to 
preference  propagation.  In  an  iterative  approach  one  cari  use  the  resource  demand  densities 
obtained  by  the  previous  iteration  to  estimate  the  probability  that  a  given  resource  will  be 
available  for  an  activity  at  some  time  r.  Using  these  probabilities,  new  a  priori  start  time 
probability  distributions  can  be  obtained  and  the  propagation  process  can  be  carried  out  all  over 
again.  Alternatively,  in  a  hierarchical  scheme,  one  can  use  the  propagation  results  obtained  at  an 
upper  level  to  compute  resource  availability  estimates.  Again  tliese  estimates  can  be  combined 
with  the  start  time  utility  functions  to  obtain  a  priori  start  time  probability  distributions  for  the 
new  level. 
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4.  Propagatiog  Start  Time  and  Duration  Distributions  in  a  TCG 


4.1.  Preliminary  Remarks 

Now  that  we  have  some  a  priori  distributions  for  the  start  time,  duration  (if  variable)  and 
resources  of  an  activity,  we  can  refine  these  a  priori  probabilities  so  as  to  account  for  the  actual 
constraints  of  the  problem.  In  this  section  we  compute  a  posteriori  start  time  (and  duration) 
probabilities  that  account  for  the  interactions  defined  by  the  TCG. 

This  section  is  subdivided  into  several  subsections  each  dealing  with  the  propagation  problem 
under  increasingly  more  general  assumptions.  As  we  already  mentioned  earlier  we  view  TCGs  as 
undirected  graphs.  Subsection  4.2  develops  the  computation  of  a  posteriori  start  time  probability 
distributions  in  acyclic  TCGs  with  rixed  duration  activities.  In  subsection  4.3  the  fixed-duration 
assumption  is  relaxed.  In  subsection  4.4  we  relax  the  acyclicity  assumption.  Subsections  4.2  to 
4.4  all  assume  that  there  are  no  explicit  disjunctions  in  the  TCG.  The  transitivity 
properties  [Alien  83]  of  the  relations  may  however  induce  disjunctions,  which  are  then  implicitly 
accounted  for  in  the  propagation.  An  exan^le  of  TCG  with  no  explicit  disjunction  is  represented 
in  Figure  4-1.  The  TCG  specifies  that  I,  OVERLAPS  Ij,  and  Ij  is  DURING  I3.  This  implicitly 
induces  the  disjunction  I,  {DURING,  STARTS,  OVERLAPS}  I3,  i.e.  Ij  is  DURING  or  STARTS 
or  OVERLAPS  I3.  All  our  calculations  allow  for  this  type  of  implicit  disjunctions.  On  the  other 
hand  explicit  disjunctioris  are  more  difficult  to  handle  and  are  quite  infrequent  in  praaical 
factory  scheduling  problems.  Propagation  in  TCGs  with  explicit  disjunctions  is  discussed  in 
subseaion  4.5.  Figure  4-2  displays  a  TCG  with  explicit  disjunctions.  Subsection  4.6  interprets 
the  results. 


explicit  constraint 
implicit  constraint 


C,;  I,  OVERLAPS 
C^:  [^DURING 

C3:  I, {DURING,  STARTS,  OVERLAPS) 


Figure  4-1:  Example  of  a  TCG  with  no  explicit  disjunctions 
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explicit  constraint 


C  :  I, {BEFORE,  AFTER) 


Figure  4-2:  Example  of  a  TCG  with  explicit  disjunctions 


Notations 

•  will  denote  the  explicit  temporal  relation  constraints  that  define  the 
TCG  (See  Figure  4-1  for  an  example). 

•  aijist^=t)  will  denote  I^’s  a  priori  start  time  probability  density,  obtained  as  suggested 
in  the  previous  section,  r  is  the  variable. 

•  wUl  denote  I^’s  a  priori  duration  probability  density,  obtained  as  suggested 
in  the  previous  section,  d  is  the  variable.  This  distribution  will  only  be  used  for 
variable-duration  activities. 

•  P(sfj.=r&C,&C2&...&C^)  (where  is  assumed  to  be  a  fixed-duration  activity) 
will  denote  I^’s  a  posteriori  start  time  probability  density.  This  density  (with 
variable  t)  corresponds  to  the  a  posteriori  probability  that  Aj.  starts  at  time  t  (i.e. 
jr^=t)  and  that  the  temporal  relation  constraints  Cj,  C2,.-,C^  are  satisfied. 

•  (where  Aj.  is  assumed  to  be  a  variable- 
duration  activity)  will  denote  the  two-dimensional  joint  a  posteriori  probability 
density  of  Ij.’s  start  time  and  duration.  This  density  (with  variables  r  and  d) 
corresponds  to  the  a  posteriori  probability  that  Aj.  starts  at  time  r  (i.e.  sx^=t)  and  has 
duration  d  (i.e.  du^=d)  and  that  the  temporal  relation  constraints  C,,  Cj, ...,  C„  are  all 
satisfied. 

In  order  to  avoid  the  accumulation  of  parentheses  in  iterated  integrals,  we  adopt  the  usual 
convention  that: 

denotes  f*(/(^)f‘^%(Ti)dTi]dq 

Ja  Ja  •'tty 

We  will  also  be  using  the  following  functions: 

•  <x(predicate)  is  a  function  that  returns  1  when  predicate  evaluates  to  true  and  0 
otherwise^. 


^Thc  reader  who  is  not  familiar  with  this  formalism  can  think  of  it  as  a  convenient  way  of  c.xprcssinc  IF- 
staicincnts  in  inailieinuiical  formulas. 
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•  P‘(EQ,(q)],  where  EQ,(^)  is  a  linear  equation  in  is  a  distribution^  such  that: 

=  a[L  <  j:  <  U\g{x)  if  x  is  the  unique  solution  to 
=j^g(^)d^  if  EQ^(^)  holds  for  V^,  and 
=  0  if  is  inconsistent. 

This  simply  expresses  that  the  integration  variable  4  is  not  only  restricted  to  values 
between  L  and  U  but  that  the  values  it  can  take  should  also  satisfy  the  linear 
equation 


'  More  generaUy  where  EQi%^ . ^„)(for  i=l  to  /)  is 

a  linear  equation,  is  a  distribution  such  that: 

. . 

glF'i^ 

if  the  system  of  linear  equations  is  consistent  and  equivalent  to: 

4,=F'(4„4.„ . t) 

. 4,) 


4)“^*  '44+j’"’’^^)  ^ 

=0  if  the  system  of  equations  is  inconsistent. 

4.2.  Propagation  in  an  Acyclic  TCG  with  Fixed>duration  Activities 
In  this  subsection  we  express  the  a  posteriori  probability  density  Pist^  =  t&C^&C2&...&  C„) 
for  the  start  time  of  an  arbitrary  time  interval  Iq  in  terms  of  the  a  priori  start  time  distributions®. 
We  assume  fixed-duration  activities  arranged  in  an  acyclic  TCG. 

We  denote  by  /J,  ...,  the  intervals  directly  adjacent  to  in  the  TCG  (Figure  4-3). 

Cf  (1  ^  i  ^  P(,)  is  the  temporal  constraint  between  I,,  and  I®.  Each  time  interval  I®  (1  <  /  <  pg)  is 
itself  related  directly  or  indirectly  to  some  other  time  intervals  by  a  set  of  constraints  S°.  The 
sets  S®  are  disjoints  as  the  TCG  is  assumed  to  be  acyclic  (Figure  4-3). 

Since  the  TCG  may  be  disconnected,  which  is  the  case  when  there  are  several  independent 
orders  to  schedule,  we  have: 


^Our  p  distribution  is  a  variation  of  the  Dirac  distribution.  The  reader  who  is  not  familiar  with  this  formalism  can 
simply  look  at  it  as  a  convenient  way  of  expressing  constraints  on  the  values  that  an  integration  variable  can  take. 

*1q  is  an  arbitrary  element  of  1 1 , ,  Ij } . 


IS 


Figure  4-3:  An  Acyclic  TCG. 

The  time  interval  In  is  related  to  /?,  fi . /° 

by  respectively  C^,  . ■ 


and  m<n-\  (number  of  edges  in  an  n-venex  tree). 

We  will  express  the  (a  posteriori)  probability  that  Iq  starts  at  time  t  and  that  the  constraints  Cj, 
are  satisfied  in  terms  of 

•  the  a  priori  probability  that  Iq  starts  at  r,  and 

•  the  probabilities  that  each  time  interval  if  (1  ^  ^-Pq)  has  a  start  time  compatible 
with  Cf  given  that  stQ=t. 

These  latter  probabilities  can  be  expressed  in  a  similar  fashion,  thereby  resulting  in  an  inductive 
formulation  of  the  a  posteriori  start  time  probabilities.  The  inductive  formulation  process  stops 
when  all  related  time  intervals  have  been  accounted  for  (or  more  precisely  their  a  priori  start  time 
distributions).  At  that  point  we  have  an  expression  of  the  a  posteriori  stan  time  distribution  that 
only  contains  a  pnori  start  time  distributions,  i.e.  distributions  that  we  know  from  the  previous 
section. 

Indeed,  the  a  posteriori  probability  that  stQ=t  and  that  die  temporal  relation  constraints  Cj. 


Ct . are  satisfied  is  given  by^  the  a  priori  probability  that  Iq  starts  at  time  f,  denoted 

GQ(stQ=t),  multiplied  by  the  conditional  probability  that  C[,  C2,...,  are  satisfied,  given  that 
stQ=t,  denoted  P  (C,  &  C,  & . . .  &  C„  =  r)): 

F(stQ  =  c&C^&C2&--. .&■€„)=  X  P(C^&.C2&.  ..&C„\srQ  =  r)  (2) 

with  r  being  the  distribution  variable. 

Furthermore,  assuming  Iq’s  start  time  fixed  at  time  r,  the  satisfaction  of  the  constraints 

{C°}u5f  is  independent  of  the  satisfaction  of  the  constraints  {C®}u5°  (for  i  ^  j),  since  we  are 
dealing  with  an  acyclic  TCG.  Hence: 

Po 

P(C^&C2&...&CJsTo  =  t)=  ^P«^&S°\st^  =  t)  (3^ 

where  P(C°&^\sIq  =  r)  is  the  conditional  probability  distribution  that  Cf  and  the  constraints  in 
S°  are  satisfied  given  that  3X^  =  1  (with  t  being  the  distribution  variable). 

Using  (3)  we  can  now  account  separately  for  each  constraint  C°.  We  will  express  each 
multiplicand,  P(C^ &.^\stQ  =  t),  in  terms  of  F(rr°  =  t&Sf),  the  probability  that  starts  at  some 
time  X  (to  be  defined)  and  that  the  constraints  in  S°  are  satisfied.  Consequently  equations  (2)  and 
(3)  win  enable  us  to  express  P(stQ  =  t&C^&C2& ■■•&€„)  in  terms  of  probabilities  of  the  form 
P(st°  =  x  &S°),  i.e.  probabilities  of  the  same  form  as  the  original  probability 

P {stQ  =  t &.C ^&.C2&. . . . &.C ^  except  that  S°  is  only  a  subset  of  (Cj,  C2,...,Cn^}.  By  recursively 
repeating  this  process,  we  will  be  able  to  account  for  all  the  temporal  relation  constraints.  The 
recursion  process  stops  when  S°  gets  empty  since  at  that  point 

F(jr°=T&5°)  =  Pistil  =  t)  =  0°(5f°  =  t),  I°’s  a  priori  start  time  density. 

Paragraghs  4.2.1  and  4.2.2  develop  the  computation  of  PiCf- t)  in  terms  of 
P(,sil  =  x&^)  =  G°(jr°=T)P(S?|sr°=t)  in  the  case  where  Cf  is  respectively  of  the  form  "Iq 
MEETS  l”",  and  "Iq  BEFORE  1°",.  The  treatment  of  the  set  of  all  thirteen  of  Allen’s  temporal 
relation  constraints  can  be  found  in  appendix  1. 

4.2.1.  C®:  Ig  MEETS  ij* 

The  constraint  '7g  MEETS  /®"  requires  (Figure  4-4)  that  /q’s  end  time  be  equal  to  /°’s  start 
time,  i.e.  it  requires  that  etQ=stQ+duQ=st°.  In  other  words,  assuming  that  stQ=t,  the  probability  that 
C°  and  the  constraints  in  S®  are  satisfied  is  equal  to  the  probability  that  st^=t+duQ  and  that  the 
constraints  in  S°  are  satisfied: 

F(C®&5®|srg  =  r)  =F(5r®  =  t  +  du^&^) 

=C5°{st°  =  t  +  dMg)  P(S®  Ur°  =  t  +  du^)  (4) 

V(/4  &  B)  =  /’(A )  X  /*(  S|A):  the  joint  probability  of  two  events  A  and  B  can  be  expressed  as  the  product  of  the 
probability  of  event  A  with  the  conditional  probability  that  B  occurs  given  that  A  is  assumed  to  occur. 
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time 


Figure  4-4;  4  MEETS 

4^.2.  Cf:  4  is  BEFORE  ij’ 


time 


Figure  4-5:  4  BEFORE 

The  constraint  "4  BEFORE  requires  (Figure  4-5)  that  /q’s  end  time  be  smaller  than  T^’s 
start  time,  i.e.  it  requires  that  etQ=stQ+duQ  <  st^ .  In  other  words,  assuming  that  s[q=[,  the 
probability  that  C°  and  the  constraints  in  Sf  arc  satisfied  is  equal  to  the  probability  that 
st°  >  t+duQ  and  that  the  constraints  in  S°  are  satisfied: 

=  =  f"  P{sP.=x&^)dx 

J  (  +  </«o 

=  f”  a°(5r”  =  T)P(5"Ur°  =  T)t/x 

•'<  *duQ 


(5) 
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4.23.  Example 

We  use  the  three  activities  of  order2  (Figure  1-2)  to  illustrate  the  computations  that  we  have 
just  developed.  We  have: 

F(srg=f&C  j  «feC2&...&C-7) 

=P{st^=t&C^&Cj) 

f  <S'jist-j=z-j)PiC‘j\st-f=X'j)dZj 

Ji+au, 

Also,  using  the  formula  given  in  appendixl  to  account  for  "I^  AFTER  Ig": 


P(,st-j=t&C  i&C2&-.-&C-j) 

=P(st-j=(&C^&C-j) 

=07(5/7=0  P(C^ScCj  I 

P(^C  P(C-j\st2~t) 

=07(5/7=0 (7) 
Finally,  in  the  same  fashion: 

/’(5/g=/&Cj&C2&.,.&C7) 

=P(st^=t&C  g&C7) 

=03(5/3=/)  F(Cg&C7 1  -^^8=0 

=03(5/3=0 1^07(5/7=^)  F(Cg  1 5/7=t7)  dij 

=03(5/3=/)  j^07(5/7=X7)dt7j^'^‘06(5/g=Xg)dt6  (8) 

We  assume  that  the  resources  are  initially  free  (i.e.  no  priori  resource  reservations).  Therefore, 
since  activities  Ag,  A7,  and  Ag  have  uniform  start  time  utility  functions,  their  a  priori  start  time 
densities  are  uniform  as  well,  and  span  between  times  0  and  120.  Figure  4-6  displays  the  a 
posteriori  start  time  densities  confuted  using  these  a  priori  densities.  Since  in  this  case  the  start 
time  utilities  are  uniform,  the  most  preferable  start  times  for  each  activity  are  the  ones  that  leave 
the  most  freedom  to  the  other  activities  for  satisfying  the  temporal  constraints  Cg  and  C7.  For 
instance,  Ag’s  a  posteriori  start  time  density  indicates  that  Ag  should  start  as  early  as  possible  in 
order  to  leave  as  much  room  as  possible  to  A7  and  Ag^®.  As  we  will  see  in  the  example  of 
subsection  4.4,  the  propagation  of  nonuniform  start  time  utility  functions  such  as  the  ones  of  the 
activities  in  orderj  are  influenced  by  a  second  factor.  In  addition  to  looking  for  start  times  that 
leave  a  lot  of  slack  to  the  other  activities  that  have  not  been  scheduled  yet,  the  propagation  of 
non'miform  utilities  gives  a  higher  preference  to  higher  utilities.  The  a  posteriori  distributions 


‘*^11  is  impoitant  at  this  point  to  bear  in  mind  that  we  have  not  yet  accounted  for  resource  capacity  constraints.  In 
particular  we  do  not  know  the  effects  that  these  constraints  will  have  on  the  domain  of  possible  start  times  for 
and  Ag.  Hence,  at  this  point,  the  best  start  times  are  the  least  committing  ones  with  respect  to  the  temporal 
constraints. 
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thereby  reflect  a  compromise  between  the  utilities  to  optimize  and  the  need  to  leave  enough 
room  for  selecting  good  stan  times  for  the  other  activities  that  have  not  yet  been  scheduled. 
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Figure  4-6:  A  posteriori  start  time  densities  for  order2 

4.3.  Propagation  in  an  Acyclic  TCG  with  Variable-duration  Activities 
In  the  case  of  variable-duration  activities,  one  can  compute  for  each  activity  a  two-dimensional 
joint  a  posteriori  probability  density  of  the  activity’s  start  time  and  duration.  This  density,  of  the 
form  PistQ=t&.duQ=d&CiScC2&—-&C„),  represents  the  probability  that  stQ=t  and  duQ=d  and 
that  all  the  temporal  relation  constraints  are  satisfied  given  the  activities’  a  priori  start  time  and 
duration  distributions.  The  integrals  involved  in  the  computation  of  these  distributions  are  very 
similar  to  those  of  the  previous  subsection,  except  that  we  now  have  to  account  for  the  a  priori 
duration  distributions. 
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The  equivalent  to  equations  (2)  and  (3)  are; 

P(,stQ  =  tScduQ=ci&Ci&C2&----&C^)=  a(stQ  =  t)  xd^idUfj^d)  x 

P{C^&C2&^^■&■C„\stQ  =  t8uiuQ=d)  (9) 

with: 

Po 

PiCi&C2&...&.CJstQ  =  t&duQ=d)=  J][/>(C:f&S°|5ro  =  r&dMo=c/)  (10) 

is  1 

Paragraghs  4.3.1  and  4.3.2  develop  the  computation  of  P(Cf  &5?| jTq  =  t&duQ=d)  in  terms  of 
P(jr°  =  t&du°  =  5&5®)  =  c.(5rf=T)5?(<iwf=5)P(5°|5r°=T&du°=6)  in  the  case  where  Cf  is 
respectively  of  the  form  "Iq  MEETS  I?",  and  "Iq  BEFORE  ij*",.  The  treatment  of  the  set  of  all 
thirteen  of  Allen’s  temporal  relation  constraints  can  be  found  in  appendixl.  Exactly  like  in  the 
previous  subsection,  one  can  use  these  equations  in  a  recursive  fashion  to  express  the  a  posteriori 
Stan  time  and  duration  densities  in  terms  of  the  a  priori  ones. 

4J.1.  Cf :  lo  MEETS 

As  before,  the  constraint  "Iq  MEETS  requites  (Figure  4-4)  that  /q’s  end  time  be  equal  to  7°’s 
Stan  time,  i.e.  it  requires  that  etQ=stQ+duQ=st°i.  Hence,  assuming  that  stQ=t  and  ,duQ=d,  the 
probability  that  C°  and  the  constraints  in  S?  are  satisfied  is  equal  to  the  probability  that  st^=t+d 
and  that  the  constraints  in  S°  ate  satisfied: 

P(C°&^lstQ  =  t&duQ=d)  s  P(sT'l=t+d&^ 

=  f  “Sf  (duf = 5)  a°(nf = f + d)  Z’  (S®  I  = r + &duf = 5)  d5 

Jo 

=jJ’5°(du°=5)d5j~3K'c=f+d)o°(jr?=x)F(S°|n°^&dM°=5)dT  (11) 

The  first  equality  is  the  most  natural  one.  We  will  however  use  equation  (11)  in  the  next 
subsection,  when  allowing  for  cycles  in  the  TCG. 

4J.2.  Cf :  Iq  before 

In  the  same  fashion,  if  Iq  has  to  be  BEFORE  I®,  one  has: 

PiC°&f^\stQ  =  t&duQ=d) 

=  f"5®(di<®=5)d5f"  <j'j(st°=x)P(^\st°=x&du^=5)dx  (12) 

Jo  J/+</ 

4.4.  Relaxing  the  Acyclicity  Assumption 

We  now  turn  our  attention  to  the  case  where  there  may  be  cycles  in  the  TCG.  Equation  (9)  still 
holds  but  the  computation  of  P(C,  =  t&dup=d)  becomes  more  complex. 

P(Ci&C2&...&C^\srQ  =  c&duQ=(f)  is  the  probability  that  C,,  Cj,  ...,  are  satisfied  when 
sto=t  and  duo=d  given  the  activities’  a  priori  start  time  and  duration  distributions.  The  set  of 


temporal  relation  constraints  {C,,  Cj,  C^}  can  be  expressed  as  a  set  of  linear  equalities  and 
inequalities  (e.g.  I,  BEFORE  is  equivalent  to  st  +du-  <  stp.  This  set  of  equalities  and 
inequalities  together  with  the  conditions  stQ=t  and  du^^d  defines  a  polyhedron  in  the  2(n-l) 
dimensional  space  generated  by  stp  st^,  st*.p  duj,  duj,  du*,,^^  The  volume  contained  in 
this  polyhedron  is  the  domain  of  admissible  values  for  st*,  st^,  st’  p  du*,  duj,  du^,  given 
the  TCG  and  the  conditions  stQ=t  and  duQ=d  (independently  of  the  a  priori  start  time  and  duration 
distributions).  Therefore  P(C^&.C2&■^^^&■C„\stQ  =  tSuiuQ=d)  can  be  obtained  by  integrating  the 
multivariable  probability  density 

aJ(st*=T,)a2(st^=t2)...CT^,(st’.,=T^,)5*(duJ=5,)62(dt4=52)...5*.,(du*.,=5^,)  over  this  volume.  In 
this  subsection  we  explain  how  to  effectively  build  this  multiple  integral  as  an  iterated  integral. 

Notice  that,  according  to  Fubini’s  theorem  (see  (Thomas  83]  for  instance),  there  axe  [2(n-l)]! 
corrert  ways  to  express  a  2(n-l)-tuple  integral  as  an  iterated  integral  (each  corresponding  to  a 
permutation  of  the  2(n-l)  integration  variables).  The  algorithm  that  we  present  builds  one  of 
these  [2(n-l)]!  iterated  integrals.  Although  all  the  iterated  forms  are  theoretically  equivalent, 
some  result  in  faster  numerical  evaluation  than  others  The  algorithm  that  we  describe  gives 
one  way  to  build  these  integrals.  The  integrals  can  then  be  rearranged  in  order  to  speed  up  their 
evaluations.  We  will  not  be  concerned  here  with  these  implementation  details. 

Consider  again  the  TCG  associated  to  order2  in  the  example  of  subsection  1.3  (Figure  4-7)^^. 

As  we  saw  in  equation  (6); 

•^(Cg&C?  I  sfg = r)  =  I"  cj2(i’r2=X2)«iT2  f  Og(srg=Tj)dTg  (13) 

Alternatively,  we  can  start  integrating  on  stg,  which  produces: 

P(C^&C2\st^=t)=  f**Og(sr8=Tg){x(Xg-du2  >  t+du^dx^  r^~^'’(j{sty=ij)dXj 

—  f  ^”^8  f  ^^(.St■J=X2)  dX^  (14) 

where  cdx^-du^  >  t+du^)  simply  expresses  that  the  second  integral’s  upper  bound  has  to  be 
greater  than  its  lower  bound  (since  we  are  integrating  probability  densities).  Figure  4-7 
represents  the  domain  of  integration  of  both  form  (13)  and  (14).  They  are  obviously  the  same, 
which  illustrates  that  both  forms  (both  iterated  integrals)  are  equivalent.  Besides  its  illustration 
of  Fubini’s  theorem,  this  example  shows  how  to  account  for  several  constraints  at  the  same  time 
when  determining  a  variable’s  domain  of  integration:  in  (14)  the  domain  of  integration  of  st^  is 
determined  by  the  two  constraints  Cg  and  C,.  In  this  example  Cg  determines  the  lower  bound  of 
the  integral  and  C^  the  upper  bound.  The  simplicity  of  the  formulas  in  the  previous  subsections 
was  coming  from  the  fact  that  it  was  possible  to  order  the  integration  variables  so  as  to  account 

“Where  {I<j,  Ij,  Ij,  ....  I„.,}  =  (Ij,  I2.  •••,  !„},  Iq  being  an  arbitrary  time  interval  of  the  set.  as  in  the  previous 
subseaions.  The  simply  indicates  that  the  time  intervals  have  been  reordered. 

'^Some  iterated  forms  are  also  easier  to  solve  analytically  than  others. 

'^One  of  the  reasons  for  chosing  this  example  is  that  the  domain  of  integration  can  be  visualized  in  2-D. 
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Figure  4-7:  Illustration  of  Fubini’s  Theorem  in  a  TCG  with  3  Time  Periods 

for  only  one  constraint  at  a  time.  In  TCGs  with  cycles  it  is  generally  not  possible  to  find  such  an 
ordering.  The  domain  of  integration  of  a  variable  is  usually  determined  by  several  constraints, 
some  affecting  the  lower-bound  some  the  upper-bound.  The  actual  lower-bound  will  therefore  be 
given  by  the  maximum  of  the  lower-bounds  produced  by  each  constraint  (i.e.  the  most  restrictive 
one)  and  the  actual  upper-bound  by  the  minimum  of  the  upper-bounds  produced  by  each 
constraint.  Additionally  one  has  to  ensure  that  the  lower-bound  is  smaller  than  the  upper-bound 
(see  a  function  in  (14))  since  we  are  integrating  probability  densities. 

Ordetj  (Figure  1-2)  in  subsection  1.3,  is  an  example  of  a  TCG  with  cycle  where  the  integration 
bounds  of  some  start  times  are  obtained  by  taking  the  minimum  or  maximum  of  the  bounds 
produced  by  several  constraints.  For  instance: 

F(C^&C2...&Cs|sr^=r)=^  aj(jr2=t2)<fT2  f 

f  a3(sr3=-C3)dt3  [  (5^{st^=x^)dx^  (15) 


As  suggested  by  the  above  exanipie,  tite  procedure  for  building  iterated  integrals  to  compute  a 
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posteriori  probability  distributions  in  TCG  with  cycles  is  just  a  generalization  of  the  formulas 
given  in  tiie  previous  subsections.  The  main  differences  come  from  the  fact  that  it  is  not  possible 
anymore  to  find  an  ordering  of  the  integration  variables  that  would  allow  for  accounting  for  only 
one  temporal  relation  constraint  at  a  time.  We  just  saw  how  to  combine  the  lower-bounds  and 
upper-bounds  imposed  by  different  temporal  constraints.  Before  describing  a  general  procedure 
to  effectively  build  a  posteriori  probability  integrals,  we  still  have  one  detail  to  consider.  Some  of 
the  formulas  given  in  the  previous  subseaion  include  3  distributions  (see  also  appendix  1).  For 
instance  p‘(T=t-Ki)  in  (1 1)  expresses  that  st?  should  be  equal  to  t-Hl  in  order  for  Iq  to  meet  I®,  p 
distributions  provide  an  easy  way  to  formally  handle  all  temporal  relation  constraints  in  the  same 
fashion,  Le.  with  integrals  over  the  duration  and  start  time  of  each  time  interval.  3  distributions 
allow  for  expressing  equalities  involving  integration  variables.  When  several  constraints 
involving  3  distributions  affect  the  same  time  interval,  one  has  to  make  sure  that  the  values  for 
the  interval’s  start  time  (or  duration)  that  they  each  require  are  compatible.  This  is  accomplished 
by  using  the  following  rule  (which  can  easily  be  verified  using  the  definition  of  3  distributions): 

J  Domain 

=  f  . (16) 

J  Domain 

This  is  illustrated  by  the  example  below. 


C,:  r,STARTED-BY  1^ 
q;  I,  CONTAINS  I, 
C^:  I  ^  MEETS 


Z  1  3 

_ ^ 

I. 

- ► 

time 

Figure  4-8:  A  TCG  with  3  Time  Periods 


The  TCG  repre.sented  in  Figure  4-8  involves  3  time  intervals  (writh  variable  durations),  namely 
I|.  Ij,  and  Ij.  Tlie  tempor;il  relation  constraints  are: 
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•  C,:  I,  STARTED-BY  I2, 

•  Cj:  I,  CONTAINS  I3,  and 

•  Cjili  MEETS  I3. 

Using  equations  (39),  (41),  (35)  (see  appendix),  and  (16)  one  can  write; 

P(C  i&Cj&Cj  I  st^=tSuiu  i=d) =j%(,du^=z^)dz^  j 

hji,dWf=^)dZ^  p2(X2=f,T2=T3-e2)a2(jt2--t2)^^'t2 

^  Max\Ofl)  Jitocl—.— ) 

This  formula  can  be  simplified  using  the  definition  of  P  distributions: 

/’(Ci&Cj&Cj  I  jrj=t&tiu,=<i) 

= J^‘'53(t/U3=e3)ife3  |^'■^“*3oJ(J^J=x3)oc(0<  tj-t  <ti)  h.^{du^=^^-t)  <52(^2=^)  ^'^3 

=  03(5X2=^)  ^{du^2~^)  dx^ 


To  conclude  this  subsection.  Figure  4-9  gives  a  description  of  BUILD-A-POSTERIORI- 
PROB  ABILITY -EXPRESSION,  a  general  procedure  to  effeaively  express 
PiC^,...,C„\stQ=tScduQ==d)  as  an  iterated  integral.  The  body  of  the  procedure  makes  use  of  a 
couple  of  simple  functions,  of  which  we  only  give  an  infoimal  description: 

•  adjacent(/,  TCG):  returns  a  list  containing  the  time  intervals  adjacent  to  /  in  the 
TCG,  I  being  itself  a  time  interval. 

•  pop(/irr):  removes  the  first  element  from  list  and  returns  it 

•  index(/):  returns  the  index  (i.e.  subscript)  of  /,  where  /  is  a  time  interval  (e.g. 
index(/2)  returns  2 ). 

•  intersection(/irr,  ,7/5X2):  returns  a  list  containing  the  elements  of  list^  that  are  also  in 
//rXj  (the  order  of  the  elements  in  the  result  list  is  arbitrary). 

•  uaion(/isti,/lst2):  returns  a  list  containing  any  element  that  is  either  in  listj  or  7/5X2  (or 
in  both).  An  element  appearing  in  both  list,  and  7/5X2  is  returned  only  once. 

•  Iist-difference(7/5X,  ,7/5X2):  returns  a  list  with  the  elements  of  7/5X,  that  are  not  in  7/5X2- 

•  start-time-upper-bound-expression(/,  list):  /  is  a  time  interval,  and  list  is  a  list  of 
time  intervals  adjacent  to  /  in  the  TCG.  The  funaion  returns  the  start  time  upper- 
bound  expression  resulting  from  the  temporal  relation  constraints  between  /  and  the 
time  intervals  in  list.  As  explained  earlier  in  this  subsection,  this  is  expressed  as  the 
minimum  of  the  upper-bound  produced  by  each  constraint. 

•  start-time-lower-bound-€xpression(/,  list):  same  as  above.  The  lower-bound  is 
expressed  as  the  maximum  of  the  lower-bound  produced  by  each  constraint. 
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•  duration-upper-bound-expression(/,  list):  same  as  above  for  the  duration 
(minimum). 

•  duration-iower-bound-expression(/,  list):  same  as  above  for  the  duration 
(maximum). 

•  beta-expression(/,  list):  combines  the  eventual  P  distributions  resulting  from  the 
temporal  constraints  between  I  and  the  time  intervals  in  list,  as  explained  earlier  in 
this  subsection.  If  there  are  no  P  distributions  the  function  simply  returns  the  empty 
expression. 

•  append(e;tpre.srion,,  expression-^:  appends  the  two  expressions  together. 

We  use  e  and  nil  to  respectively  denote  the  empty  expression  and  the  empty  list. 


procedure  BUILD-A-POSTERIORI-PROBABIUTY-EXPRESSION  (Tq  ,  TCG) 

INTERVALS-TO-BE-PROCESSED  <-  adjacent(Io,  TCG) ; 

PARTIAL-EXPRESSION  £ ; 

MARKED-INTERVALS  <-  {Iq}  ; 

*a  list  containing  Ig* 

while  INTERVALS-TO-BE-PROCESSED  ^  nil 

I  <-  pop(INTERVALS-TO-BE-PROCESSED)  ; 
i  <-  index(I) ; 

RELATED-INTERVALS  adjacent(I,  TCG) ; 

INTERVALS-TO-ACCOUNT-FOR  «- 

intersection(MARKED-INTERVALS,  RELATED-INTERVALS)  ; 
SUB  start<time<upper>bound<expression(I,  iNTERVALS-TO-ACCOUNT-FOR) ; 

SLB  4-  start-time-lower<bound>expressionG,  INTERVALS-TO-ACCOUNT-FOR)  ; 

DUB  4—  duration«upper>bound-expression(I,  INTERVALS-TO-account-FOR)  ; 

DLB  4—  duration-lower-bound>expression(i,  INTER valS-TO-ACCOUNT-FOR)  ; 

BETA  4-  beta-expression(I,  intervals-TO-ACCOUNT-FOR)  ; 

LOCAL-EXPR  4-  a(DC/B>DLS)J“'"5,<du=e,)<X(SUR>SLR)de,.|*'"fl£rAa,(jr..=T,.)dT,. ; 

PARTIAL-EXPRESSION  4-  append(PARTIAL-EXPRESSIONJLOCAL-EXPR)  ; 
MARKED-INTERVALS  4-  union(MARKED-INTERVALS,  {I})  ; 
INTERVALS-TO-BE-PROCESSED  4-  union(INTERVALS-TO-BE-PROCESSED, 

list-difference(RELATED-INTERVALS,  MARKED-INTERVALS))  ; 

while-end ; 

return  PARTIAI^EXPRESSION; 


Figure  4-9:  Procedure  to  express  P{Ci,...,CJstQ=tSuiuQ=d)  as  an  iterated  integral 

The  procedure  builds  the  iterated  integral  from  left  to  right  by  successively  visiting  each  time 
interval  that  is  directly  or  indirectly  related  to  Ig.  LOCAL-EXPR  contains  the  integrals  over  the 
start  time  and  duration  of  the  time  interval  currently  visited.  This  local  expression  is  appended  to 
the  right  of  a  current  partial  expression  of  the  iterated  integral,  thereby  resulting  in  a  new  partial 
expression.  Intervals  that  have  been  visited  (i.e.  whose  start  time  and  duration  a  priori  densities 
have  already  been  integrated  in  PARTIAL-EXPRESSION)  are  marked.  Integration  bounds  for  a  time 
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interval’s  stan  time  and  duration  are  determined  by  the  temporal  relation  constraints  between 
that  time  interval  and  the  adjacent  time  intervals  that  have  already  been  marked. 

As  illustrated  in  the  previous  examples,  the  expressions  produced  by  this  procedure  can  be 
simplified  using  the  definitions  of  a  functions  and  P  distributions.  The  integration  bounds  can 
also  be  refined  to  account  for  the  very  domain  over  which  the  probability  densities  are  strictly 
positive.  Finally  the  order  of  integration  can  be  rearranged  to  speed  up  evaluation.  A  time 
complexity  analysis  of  the  method  and  a  discussion  of  available  methods  to  evaluate  the  integrals 
are  given  in  section  6. 

Figure  4-10  illustrates  the  operation  of  the  procedure  in  the  construction  of  the  iterated  integral 
in  (15).  Notice  that  the  a  expressions  have  been  omitted  as  they  trivially  evaluate  to  1. 

Figure  4-11  displays  the  a  posteriori  start  time  densities  of  the  activities  in  order  j,  assuming  no 
prior  resource  reservations.  The  start  time  utility  functions  are  the  ones  described  in  subsection 
1.3,  triangle  sh^d  utility  functions  allowing  for  start  times  between  0  and  140  with  a  peak  in 
120.  One  should  notice  the  difference  with  the  propagation  of  the  uniform  start  time  utilities  of 
order2  (Figure  4-6).  For  instance  in  the  case  of  Aj,  the  a  posteriori  density  was  not  totally 
pushed  to  the  right.  Instead  the  density  peaks  at  130,  which  is  a  compromise  between  the  optimal 
start  time  (120)  and  the  tendency  of  the  other  activities  to  push  Aj  towards  its  latest  start  time 
(140)  in  order  to  have  more  freedom.  A  similar  remark  applies  to  the  other  four  activities.  It 
should  also  be  noted  that  Aj’s  a  posteriori  start  time  density  between  40  and  50  and  Aj’s 
between  90  and  100  are  not  zero,  though  very  small,  which  unfortunately  does  not  appear  very 
clearly  on  the  graphs. 
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initialization; 

PARTIAL-EXPRESSION;  E 
MARKED-INTERVALS;  {Ii} 

INTER VALS-TO-BE-PROCESSED;  {12,14} 


stepi; 

I;l2 

PARTIAL-EXPRESSION;  <52ist2=t2)^h 

MARKED-INTERVALS;  (1 1,12} 

INTER  VALS-TO-BE-PROCESSED;  {14,13} 

step2; 

I:l4 

PARTIAL-EXPRESSION;  f“  a2(st2=r2^'hl“  (y4(STi=r4)dX^ 

JfWuj 

MARKED-INTERVALS;  {Ij,  12,14} 

INTER  VALS-TO-BE-PROCESSED;  {13,15} 


step3; 

I;l3 

PARTIAL-EXPRESSION;  a2(^r2=^2)^^ f ”  <^4^^h^4)(^'^4\°° 

J(WU|  Jj-Kiuj  JtjWu2 

MARKED-INTERVALS;  {Ii,  I2, 13, 14} 

INTER  VALS-TO-BE-PROCESSED;  {I5} 


step4: 


I;l5 

PARTIAL-EXPRESSION:  a,(5r,=T,)dT,  f" 


MARKED-INTERVALS:  {I^,  I2, 13, 14, 15} 
INTERVALS-TO-BE-PROCESSED:  nil 


<y4(-yr4=T4)  d\  f  “  a3(ir3=T3)  d-Zj 

Jt,'H<U2 


f"  Z5i{st^=x^)dx^ 

JMaxi-CjHiuyX^'Hiu^] 


Figure  4-10:  Main  steps  involved  in  the  construction  of  (15). 

Notice  that  the  a  expressions  have  been  omitted  as  they  trivially 

evaluate  to  1. 
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Figure  4-11:  A  posteriori  start  time  densities  for  order  j 
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4.5.  Propagation  in  a  TCG  with  Explicit  Disjunctions 

In  the  case  of  explicit  disjunctions  in  the  TCG,  one  has  to  add  the  probabilities  of  all  possible 
combinations  of  relations.  Consider  the  example  displayed  in  Figure  4- 1 2.  There  are  three  time 
intervals  I,,  I^,  and  I3.  The  temporal  relation  constraints  are; 

•  C,:  I,  (BEFORE,  AFTER}  Ij,  and 

•  Cj;  I2  (BEFORE,  AFTER}  I3. 

Four  combinations  are  possible: 

1. 1,  BEFORE  Ij  and  BEFORE  I3, 

2. 1,  BEFORE  I2  and  AFTER  I3. 

3.  Ij  AFTER  Ij  and  Ij  BEFORE  I3,  and 

4. 1,  AFTER  and  AFTER  I3. 

The  a  posteriori  probability  distributions  are  obtained  by  adding  the  probabilities  of  all  four 
possibilities.  For  example; 

P(Ci&C2lst,=t) 

Jr+au| 

+  f  ‘~^-ky^(st^=z-^)dX2  f  **  a3(5r3=T3)dT3 

J  J  t2+au2 


C,:  I,  {BEFORE,AFTER} 
q:  (BEFORE, AFTER}  I3 

Figure  4-12:  A  TCG  with  explicit  disjunctions 

As  the  number  of  possible  combinations  grows  exponentially  with  the  number  of  time 
intervals,  the  computations  are  expected  to  quickly  become  intractable.  Fortunately,  in  the 
factory  scheduling  domain,  it  has  been  our  experience  that  such  disjunctions  are  extremely 
infrequent. 
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4.6.  Result  Interpretation 

It  is  interesting  at  this  point  to  look  back  at  the  desiderata  and  requirements  identified  in 
subsection  3.2  and  check  if  they  are  satisfied.  From  subsection  3.2.  we  already  know  that  our  a 
priori  start  time  distributions  have  been  built  to  satisfy  requirement!.  In  order  to  properly 
perform  consistency  checking  we  still  need  to  check  requirement!,  i.e.  we  need  to  make  sure 
that  no  admissible  value  may  receive  a  zero  a  posteriori  probability.  This  just  follows  from 
probability  theory.  Moreover  the  method  will  give  a  zero  a  posteriori  probability  to  any  value 
forbidden  by  the  TCG,  given  the  a  priori  probability  distributions.  Hence  the  computation  of  the 
a  posteriori  probabilities  is  perfect  with  respect  to  desideratum!,  as  far  as  the  interactions 
defined  by  the  TCG  are  concerned.  The  remaining  inconsistencies  result  from  the  difficulty  to 
account  for  inter-order  interactions  between  unscheduled  activities. 

Once  the  a  posteriori  start  time  and  duration  distributions  have  been  computed,  one  has  to 
distinguish  between  two  possible  situations: 

1 .  If  at  least  one  of  the  a  posteriori  probability  density  is  uniformly  zero  then  the 
current  CSP  is  unsatisfiable  (inconsistent).  The  incremental  scheduler  should 
backtrack,  if  still  possible. 

2.  Otherwise,  after  having  been  normalized,  the  a  posteriori  distributions  can  be 
combined  to  obtain  the  resource  demand  densities  induced  by  the  CSP,  as  we 
describe  in  the  next  section^'^.  The  normalization  simply  expresses  that  the  total 
probability  that  each  activity  occurs  is  equal  to  one. 

Because  they  account  for  the  interactions  defined  by  the  TCG,  a  posteriori  start  time  (and 
duration)  distributions  reflect  intra-order  interactions.  They  generalize  the  Operations  Research 
notion  of  activity  slack  [Johnson  74].  Indeed  a  value  with  a  high  a  posteriori  probability  will 
usually  correspond  to  a  high  utility  and  will  be  likely  to  leave  a  lot  of  freedom  for  selecting  high 
utility  values  for  the  other  variables  that  have  not  been  assigned  a  value  yet.  Therefore  selection 
of  start  times  (and  durations)  with  high  a  posteriori  probabilities  is  expected  to  result  in  good 
solutions  to  the  CSP  (desideratuml)^^.  An  activity  whose  range  of  admissible  start  times  (and 
durations)  with  high  a  posteriori  probabilities  is  very  wide  is  an  activity  with  a  lot  of  slack  (with 
respect  to  the  TCG).  On  the  other  hand,  if  the  range  of  admissible  values  with  high  a  posteriori 
probabilities  is  small,  the  activity  has  little  slack.  Equivalently,  from  a  constraint  satisfaction 
point  of  view,  these  a  posteriori  distributions  can  be  seen  locally  as  measures  of  value  goodness 
and  globally  as  measures  of  variable  looseness. 

However  it  is  important  to  understand  that,  in  general,  the  peaks  of  the  a  posteriori  start  time 
and  duration  distributions  will  not  exactly  coincide  with  the  optimal  activity  start  times  and 
durations  of  the  problem,  nor  will  they  even  coincide  with  those  of  the  problem  obtained  by 
omitting  the  resource  capacity  constraints.  For  instance,  in  the  case  of  order  the  optimal  start 
times  of  Aj,  A2,  A3,  A4,  and  A5  are  respectively  30,  60,  90,  60,  and  120  (and  30,  60,  90,  90,  and 
120  if  one  omits  R2’s  capacity  constraint).  Obviously  these  optimal  start  times  do  not  exactlv 
coincide  with  the  peaks  of  the  distributions  displayed  in  Figure  4-11.  TTiis  is  because  at  this  stage 


'‘‘This  situation  is  not  a  guarantee  that  the  current  CSP  is  satisfiable  since  we  have  only  performed  partial 
consistency  checking. 

'^Although  one  should  still  account  for  inter-order  interactions,  which  is  the  topic  of  the  next  .section. 
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we  have  not  accounted  precisely  for  the  interactions  induced  by  the  capacity  constraints.  Our  a 
priori  distributions  accounted  only  implicitly  for  the  existence  of  these  interactions  by  assuming 
non  zero  probabilities  for  values  that  were  not  locally  optimal  (see  Figure  3-1).  Iterating  the 
propagation  process  as  suggested  in  subsection  3.2,  i.e.  using  resource  demand  densities  to  guess 
new  a  priori  probabilities,  should  improve  the  quality  of  the  a  posteriori  start  time  and  duration 
distributions  as  measures  of  start  time  and  duration  goodness  (desideratuml). 

5.  Resource  Demand  Densities 

We  complete  the  propagation  process  by  combining  the  a  posteriori  start  time  (and  duration) 
probabilities  to  estimate  the  amount  of  contention  for  each  resource.  This  is  performed  in  two 
steps: 

1.  For  each  activity  we  compute  a  set  of  individual  demand  densities  D/^j.  For 

each  resource  that  an  activity  can  use,  the  demand  density  D/^ji^t)  reflects 
the  probability  that  Aj^  uses  at  time  t  to  fulfill  its  resource  requirement 
This  probability  depends  both  on  the  probability  that  A^  is  active  at  time  t  and  the 
probability  that  uses  R/^j  to  fulfill  its  requirement  R/^.  The  probability  that  an 
activity  is  active  at  some  time  t  is  given  by  the  probability  that  the  activity’s  start 
time  and  duration  are  such  that  the  activity  does  not  start  after  t  and  does  not  start 
so  early  that  it  is  already  finished  by  f.  In  the  case  of  a  fixed-duration  activity  A^. , 
this  is  the  probability  that  the  aaivity  starts  some  time  between  t-dUj,  and  r.  A 
detailed  treatment  of  ±e  computation  of  individual  demand  densities  is  given  in 
appendix!.  We  will  also  interpret  as  the  reliance  of  A^  on  the  possession  of 

Rj^j  at  time  t.  Indeed  activities  with  little  slack  and  few  good  possible  resources 
will  have  high  individual  demand  densities  concentrated  over  short  time  periods 
and  a  few  resources,  whereas  activities  with  a  lot  of  slack  and  sevend  good 
resource  alternatives  will  have  smoother  individual  demand  densities  spread  over 
long  periods  of  time  and  several  resources  (see  appendix  2  for  details). 

2.  For  each  resource,  activities’  individual  demand  densities  arc  combined  to  obtain 
the  resource’s  aggregate  demand  density. Tim  density  gives  the  expected  demand 
for  the  resource  as  a  funaion  of  time.  In  the  example  described  in  subsection  1.3, 
the  aggregate  demand  densities  are  given  by: 

•  Rj’s  aggregate  demand  density  =  D^(f)=Dm(r)+Z)2n(^)+^51lW 

•  R2’s  aggregate  demand  density  =  D2(t)=Dii2(^)+I53 12(0+1^4 

•  R3’s  aggregate  demand  density  =  f?3(f)=Dgj3(r)-i-£)7j3(f)+Dg23(f) 

Notice  that  the  aggregation  process  is  performed  regardless  of  the  resources’ 
capacities.  As  a  matter  of  fact,  a  resource’s  aggregate  demand  density  at  some  time 
t  may  get  larger  than  its  capacity.  In  general  high  contention  for  a  resource  will 
require  prompt  attention  from  the  scheduler. 

Figure  5-1  depicts  the  aggregate  demand  densities  Dj(t),  D2(t),  and  D3(t)  for  the  example  in 
subsection  1.3.  Clearly  the  contention  between  A3,  A4  and  A7  for  R2,  which  was  predicted  in 
the  introduction,  has  been  identified  by  the  propagation  method.  It  corresponds  to  the  peak  of 
D2(t)  centered  around  t=100.  This  peak  reaches  a  density  of  1.5,  which  is  much  larger  than  any 
of  the  other  peaks.  Figure  5-2  shows  the  individual  contributions  of  A3,  A4,  and  A7  to  the 
demand  around  the  peak,  namely  0317(0,  0417(0,  and  0717(0.  An  area  of  width  30  has  been 
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Figure  5-1:  Rj,  R2,  and  R3’s  aggregate  demand  densities 

delimited  around  the  peak^^.  This  is  the  area  of  high  contention  for  R2.  It  clearly  appears  that 
within  that  zone,  A3  is  the  activity  whose  individual  demand  density  contributes  most  to  the 
demand  for  R2.  Consequendy  A3  is  the  activity  that  relies  the  most  on  the  possession  of  R2 
within  the  area  of  high  contention.  An  incremental  scheduler  can  accordingly  decide  to  first 


'*There  is  no  particular  reason  for  chosing  30  except  that  it  seems  to  be  a  characteristic  duration  for  this  problem, 
since  all  the  activities  have  a  duration  of  30.  The  same  results  would  hold  if  we  were  considering  slightly  smaller  or 
larger  intervals  of  contention  around  the  peak. 
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Figure  5-2:  Contributions  of  A3,  A4,  and  A7  to  R2’s  aggregate  demand  density 
focus  its  attention  on  the  scheduling  of  A3. 
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6.  Discussion 
6.1.  Time  Complexity 

As  demonstrated  in  section  4,  in  a  TCG  with  no  explicit  disjunctive  constraints, 
P{st,^  =  t8cdui=d &C,  & Cj & ...  &  C„)  can  be  expressed  at  worst  with  a  2(n-l )-tuple  integral  of  a 
priori  start-time  and  duration  probability  densities,  where  n  is  the  number  of  activities  to 
schedule.  The  construction  of  these  integrals  can  be  performed  in  polynomial  time  (see  the 
procedure  in  Figure  4-9).  On  the  other  hand,  evaluation  of  multiple  integrals  using  classical 
integration  techniques  requires  exponential  time.  In  the  worst  case  computation  of  the  n  a 
posteriori  distributions  requires  0(nK^)  integrand  evaluations,  where  K  is  a  constant  that 
depends  on  the  integration  method.  In  the  case  of  fixed-duration  activities  this  complexity  is  still 
0(nK").  This  exponential  worst-case  time  complexity  is  actually  a  very  pessimistic  one.  In 
manufacturing  environments  activities  are  grouped  in  orders.  Only  activities  within  the  same 
order  have  temporal  relation  constraints  between  them.  Therefore  the  largest  multiple  integrals 
that  one  has  to  evaluate  correspond  to  the  largest  number  of  interconnected  activities  within  an 
order  (say  Max  This  results  in  a  worst  case  time  complexity  of  nK?^“”order  in  the  case  of 

n  variable-duration  activities,  and  a  time  complexity  of  nK^”order  in  the  case  of  n  fixed- 
duration  activities.  This  also  means  that,  for  a  given  set  of  order  types  (i.e.  Max  is  fixed), 
the  worst-case  time  complexity  to  compute  the  a  posteriori  probability  distributions  is  linear  in 
the  number  of  orders  to  schedule.  The  computation  of  the  resource  demand  densities  requires  at 
most  0{n  X  r)  steps,  for  n  activities  and  r  resources.  Hence  for  a  fixed  set  of  order  types  and  a 
fixed  set  of  resources,  the  asymptotic  time  complexity  of  the  approach  is  linear  in  the  number  of 
orders  to  schedule. 

In  manufacturing  environments,  one  may  have  to  schedule  up  to  several  thousands  of  activities 
grouped  in  orders  of  up  to  20  or  30  activities.  Assuming  that  half  of  these  activities  are  modeled 
as  variable-duration  activities,  one  may  have  multiple  integrals  of  dimension  up  to  60.  Numeric 
evaluation  of  such  integrals  is  usually  performed  using  Monte  Carlo  techniques  [Stroud  71]. 
[Lepage  78]  describes  an  adaptive  Monte  Carlo  method  for  evaluating  multidimensional 
integrals  whose  asymptotic  time  and  space  complexities  are  linear  in  the  integral’s  dimension. 

Alternatively  one  may  try  to  reduce  the  size  of  the  integrals  via  the  use  of  a  hierarchical 
scheduler. 


6.2.  Expressiveness  of  the  Model 

The  preference  propagation  techniques  that  we  have  presented  allow  for  all  thirteen  of  Allen’s 
temporal  relation  constraints  as  well  as  for  disjunaions  of  such  constraints.  Additionally 
quantitative  temporal  relation  constraints  such  as  ” Activity g  should  start  at  least  5  minutes  after 
activity^"  can  be  represented  using  dummy  activities.  For  instance,  one  can  introduce  a  dummy 
activity^  with  duration  of  5  minutes  and  the  two  constraints  " activity MEETS  activity q”  and 
"activity^  BEFORE  activity g  \  Using  duration  preferences  one  can  express  even  more  complex 
quantitative  temporal  relation  constraints  such  as  "Activity^  should  start  as  soon  as  possible 
within  5  minutes  after  activity^". 

Our  model  accounts  for  three  types  of  local  preferential  constraints:  start  time,  duration,  and 
resource  preferential  constraints.  End  time  preferential  constraints  can  be  expressed  using 
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dummy  activities.  For  instance  an  end  time  preferential  constraint  on  an  activiry^  can  be 
expressed  as  a  start  time  preferential  constraint  on  a  dummy  activity^  MET-BY  activity^.  Our 
framework  also  seems  to  allow  for  the  representation  of  the  most  common  global  organizational 
constraints  [Baker  74].  For  instance,  minimizing  mean  (weighted)  order  tardiness  can  be 
expressed  with  the  help  of  end  time  (hence  stan  time)  constraints  on  the  last  activities  of  each 
order.  Minimization  of  mean  (weighted)  order  flowtime  can  be  represented  with  aggregate 
activities,  each  containing  ail  the  activities  in  an  order  and  a  preferential  constraint  on  the 
duration  of  each  aggregate  activity. 

63.  Possible  Improvements 

The  preference  propagation  technique  presented  in  this  paper  has  been  implemented  on  a  Sun 
3/60  running  Knowledge  Craft  on  top  of  Lucid  Common  Lisp  for  TCGs  with  fixed-duration 
activities  intercormeaed  by  BEFORE/AFTER  relations  and  that  may  contain  cycles.  An 
incremental  scheduler  has  also  been  built  that  uses  the  preference  propagation  module  to  focus 
its  attention.  Preliminary  experimentation  with  the  system  suggests  several  possible 
improvements. 

6  J.l.  Iterative  and  Hierarchical  Preference  Propagation 

As  already  mentioned  in  subsection  3.2,  we  are  currently  investigating  alternative  ways  to 
compute  a  priori  probability  distributions.  In  particular  we  are  considering  both  iterative  and 
hierarchical  variations  of  the  preference  propagation  scheme  presented  in  this  paper.  In  an 
iterative  approach  one  can  use  the  resource  demand  densities  obtained  by  the  previous  iteration 
to  estimate  the  probability  that  a  given  resource  will  be  available  for  an  activity  at  some  time  r. 
These  probabilities  can  then  be  combined  to  obtain  more  accurate  start  time,  duration,  and 
resource  a  priori  probability  distributions  for  a  new  propagation.  For  instance  good  start  times  for 
which  good  resources  ate  likely  to  be  unavailable  would  see  their  a  priori  probability  being 
reduced.  A  hierarchical  tq}proach  is  similar  except  that  the  additional  information  is  obtained 
from  the  results  of  the  propagation  at  the  upper  level  rather  than  from  the  previous  iteration. 
Such  techniques  are  expeaed  to  account  mote  accurately  for  the  resource  requirement 
interactions  of  unscheduled  activities. 

6  J.2.  Activity  Criticality 

In  the  introduction  we  have  defined  a  critical  activity  as  one  whose  good  (overall)  start  times 
and  resources  are  likely  to  become  unavailable  if  one  started  scheduling  other  activities  first.  In 
this  paper  we  have  assumed  that  the  most  critical  activity  is  the  one  that  relied  the  most  on  the 
possession  of  the  most  contended  resource  (over  the  area  of  high  contention  for  that  resource). 
This  measure  of  activity  criticality  is  only  concerned  with  the  availability  of  good  resources  at 
good  start  times.  Good  start  times  may  however  become  unavailable  just  because  of  operation 
precedence  interactions  (i.e.  intra-order  interactions),  as  reflected  in  the  a  posteriori  start 
time/duration  distributions.  ^^We  are  looking  for  ways  to  integrate  the  notion  of  start 
time/duration  looseness  identified  in  subsection  4.6  directly  into  our  measure  of  activity 


*^In  the  approach  that  we  have  presented,  incra-otder  interactions  are  accounted  for  indirectly  via  the  individual 
and  aggregate  demand  densities,  since  these  densities  are  computed  from  the  a  posteriori  start  time/duration 
distributions. 


criticality  rather  than  only  indirectly  through  measures  of  resource  contention  and  activity 
resource  reliance.  Additionally,  rather  than  simply  accounting  for  activity  reliance  with  respect 
to  the  most  contended  resource,  we  would  like  to  develop  a  measure  that  accounts  for  the 
reliance  of  an  activity  on  each  of  its  possible  resources  and  the  contention  on  each  of  these 
resources  (over  the  appropriate  time  intervals). 

633.  Value  Goodness 

All  along  we  have  assumed  that  value  goodness  was  solely  determined  by  the  problem 
constraints,  i.e.  both  the  required  and  preferential  constraints  of  the  problem.  A  more 
sophisticated  approach  would  consist  in  also  accounting  for  the  time  available  to  come  up  with  a 
schedule.  If  there  is  very  litde  time  available,  one  will  be  mainly  concerned  with  finding  an 
admissible  schedule  as  soon  as  possible.  Good  values  are  therefore  the  ones  that  are  the  least 
likely  to  result  in  backtracking,  i.e.  the  least  constraining  values  identified  in  earlier  work  in 
constraint  satisfaction  with  uniformly  preferred  values  [Haralick  80].  Instead  if  more  time  is 
available,  it  may  be  worthwhile  considering  riskier  values  because  they  are  likely  to  result  in  a 
better  schedule.  For  instance,  if  one  machine  is  more  accurate  than  all  the  other  ones,  one  could 
try  to  schedule  more  activities  on  the  most  accurate  machine.  This  may  however  result  into  some 
extra  backtracking  due  to  the  higher  contention  for  the  accurate  machine. 

7.  Summary  and  Concluding  Remarks 

Factory  scheduling  is  subject  to  a  wide  variety  of  preferential  constraints  such  as  meeting  due 
dates,  reducing  order  flowtime,  using  accurate  machines,  etc.  These  local  a  priori  preferences 
interaa.  For  instance,  meeting  an  order’s  due  date  may  prevent  the  scheduler  from  selecting  an 
accurate  machine  for  an  operation.  Therefore  selecting  start  times  or  resources  based  solely  on 
such  preferences  is  likely  to  result  in  poor  schedules.  Preference  propagation  strives  for  the 
construction  of  measures  that  reflect  preference  interactions.  Such  measures  can  then  serve  to 
guide  the  construction  of  good  overall  schedules  rather  than  schedules  that  locally  optimize  a 
subset  of  a  priori  preferences. 

Our  approach  to  preference  propagation  is  inspired  by  two  CSP  look-ahead  techniques  known 
as  variable  ordering  and  value  ordering  [Dechter  88].  Both  theoretical  and  empirical  studies 
[Haralick  80,  Freuder  82,  Nudel  83,  Purdom  83,  Stone  86]indicate  that  these  techniques  can 
significantly  reduce  the  amount  of  search  for  a  solution.  Earlier  work  had  only  focused  on 
applying  these  techniques  to  CSPs  where  variables  have  finite  sets  of  equally  preferred  values. 
Our  approach  to  preference  propagation  extends  these  techniques  to  CSPs  where  variables  have 
infinite  bounded  sets  of  possible  values  with  non-uniform  preferences.  The  resulta  of  the 
propagation  are  formulated  as  a  set  of  texture  measures.  In  this  paper  we  have  identified  the 
following  texture  measures:  start  time/duration  goodness  and  looseness,  resource  contention, 
and  activity  resource  reliance. 

From  an  Operations  Research  point  of  view  our  preference  propagation  technique  combines 
advantages  of  both  order-based  and  resource-based  scheduling  by  accounting  for  both  intra-order 
and  inter-order  interactions  [Smith  85]. 

We  perform  preference  propagation  within  a  probabilistic  framework.  A  probability  is 
associated  to  each  variable’s  possible  value  that  dynamically  reflects  the  likelihood  that  the  value 
results  in  a  good  schedule  overall.  We  have  identified  requirements  and  desiderata  to  guide  the 
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construction  of  such  probabilities.  These  requirements  and  desiderata  have  been  motivated  by  a 
double  objective: 

1.  We  want  to  be  able  to  detect  unsatisfiable  CSPs  as  soon  as  possible  (quick 
pruning),  and 

2.  we  want  to  use  the  propagation  results  to  help  focus  the  scheduler’s  attention  on 
the  most  critical  decision  points  and  the  most  promising  decisions  at  these  points 
(opportunistic  scheduling). 

We  have  argued  that  the  approach  presented  in  this  paper  fulfills  these  requirements  and 
desiderata. 

We  have  described  an  algorithm  to  perform  preference  propagation  in  T/CCGs.  The  algorithm 
deals  with  aU  thirteen  of  Allen’s  temporal  relation  constraints  and  allows  for  cycles  in  the 
corresponding  TCG.  The  algorithm  also  allows  for  activity  start  time,  duration,  and  resource 
preferences  and  accounts  for  earlier  resource  reservations  if  any.  We  have  shown  that  the  results 
of  the  propagation  across  the  temporal  constraints  can  be  combined  to  estimate  resource 
contention  and  activity  resource  reliance.  We  have  also  analyzed  the  computational 
requirements  of  our  approach. 

The  importance  of  this  research  lies  in  its  attempt  to  give  a  more  formal  characterization  of  the 
problem  space,  in  which  we  carry  the  search  for  a  schedule.  Given  the  underlying  uncertainty  of 
any  search  problem,  a  probabilistic  characterization  is  a  very  attractive  one.  In  this  paper,  we 
have  presented  a  model  that  uses  Bayesian  probabilities  to  account  for  preference  interactions  in 
T/CCGs.  The  problem  space  is  finally  characterized  by  a  set  of  textures  that  are  used  to  guide  the 
search  process. 
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Appendixl:  A  Posteriori  Start  Time  and  Duration  Distributions 

In  this  appendix  we  summarize  the  essential  formulas  developed  in  subsections  4.2  and  4.3  and 
complete  them  to  allow  for  all  13  of  Allen’s  temporal  relation  constraints.  The  notations  are  the 
ones  defined  in  subsection  4.1. 


1.  Acyclic  TCG  with  fixed-duration  activities 

We  found  in  subsection  4.2  that: 

P{st^  =  r&C,  &C2& . . .  &CJ=  <Jo(4ro  =  r)  x  />(C,  . .  .ScCJst^  =  t) 

with: 

pq 

„  I  . 

C“  may  be  any  of  Allen’s  thineen  temporal  relation  constraints: 

1.1.  C°:  Iq  meets  IJ* 

P(Cf&Sflsro  =  r)  =  r  +  r/«o&5f) 

=o°(sr“  =  r  +  diL^)P{f.\4  =  t  + 


1.2.  Cf:  lo  MET-BY  if 

/»(Cf&5f|sro  =  r)  =P(st°  =  t-du°,&^) 

= of  ( jr?  =  r  -  du°)  />(Sf  I  srf  =  r  -  rfuf) 

1 J.  Cf;  !<,  BEFORE  if 

P(Cf &Sf  Isfo  =  0  =  1"  F(sff  =  x&Sf)(/T 
=  r  of(stf  =  T)/>(S^jrf  =  x)t/t 

•  /  +  dUQ 


1.4.  Cf:  lo  AFTER  if 

J} 

F(Cf&Sf|5ro  =  r)  =  J‘~  ‘/’(srf  =  x&sf)<fx 
=  J'■‘^‘af(srf  =  x)/’(Sfls^f  =  x)dx 


1.5.  Cf;  lo  DURING  If 

/»(Cf  &  5f  I  sto  =  r)  =  cx(t/Uo  <  c/uf )  f  o^(srf  =  x  &  5f )  rfx 

^  f  +  du^  —  ifu  ■ 

=  aidu^  <  du^)  f '  =  X)  F(Sf  i  srf = x)  dx 

^  f  *  i/ttn  —  du : 


(17) 

(18) 


(19) 


(20) 


(21) 


(22) 


(23) 
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1.6.  C":  Iq  contains  I® 

=  aiduo  >  du°)  ^  a°(4  =  T)/’(S®|.rrf=T)^fT 

1.7.  Cf:  lo  STARTS  ij* 

/»(Cf &S°l«o  =  r)  =  a(i/uo  <  du°)  P(st°  =  r &Sf) 

=  a^duQ  <  du°)  =  f) 


1.8.  Cf:  lo  STARTED-BY  if 

/»(C®&Sf  I^To  =  f)  =  aj(duo  >  du^)  P(st°  =  r&sf) 
=  a(Ji/o  >  af(^rf  =  /)  P(Sf  l^rf =/) 


1.9.  Cf:  lo  FINISHES  if 

F(Cf  &  Sf  I  ^To  =  f )  =  (x(i/Uo  <  cfuf )  =  f  +  dtto  -  &  -^ ) 

=a(duo  <  du°)  af(5rf  =r+tiuo~J«f)f  (Sf  Urf =r  +  dup  -  du°) 

1.10.  Cf :  lo  FINISHED-BY  if 

F(Cf &Sf  liTo  =  f)  =  a(d«o  >  duj)  Pist^f  =  t  +  du^-  d«f &sf) 
= aW«o  >  d«f)  af(jff  =t+duo-du^P{f,  1 5rf = f  +  dMo  - 


1.11.  cf:  lo  OVERLAPS  if 

P(Cf&Sfl«o  =  r)  =  J 


•  dUQ 


Afox  (r,  t  +  ^un  -  du- ) 


„P(5rf  =  T&sf)rfT 


=I 


f  +  dun 


„a>f  =  T)P(Sfkrf=T)(/T 

Max  { t,  /  .)■  duQ  -  du- ) 


1.12.  cf:  lo  OVERLAPPED-BY  if 


1  -  <iu 
0  ' 

r Min  It,  I -t- dUfs  -  du ■  \  „ 

=  f  ^  '  of(5rf  =  x)P(Sf|jff=-c)^iT 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


1.13.  Cf:  Iq  equals  I® 

P(C®&5®|5ro  =  r)  =  aidiiQ=du^)Pis/!=t&^) 

= aiduQ=du^  a®(^r®=r)  P(^  \  st^=t)  (3 1 ) 


2.  Acyclic  TCG  with  variable-duration  activities 
We  found  in  subsection  4.3  that: 


P{stQ  =  t8cduQ-d&C^&C2&■-■^&■C„)=  c{stQ  =  t)  x5o(dMg=d)  x 
P{C^&C2&...&.CJstQ  =  tSuiuQ=d) 


(32) 


with: 

Po 

/>(C,&C2&...&C„Uro  =  ?&^«o=^=  YlP{C°&^\sto  =  t&dUo=d)  (33) 

1=  1 

C®  may  be  any  of  Allen’s  thineen  temporal  relation  constraints: 

2.1.  Cf:  lo  MEETS  if 

PiC°&^\st^  =  t&du^=d)  s  Pist°=t+d&^ 

=  f'‘5f(d«°=5)af(5rf=f+dO/’(Sf  Ur°=^+^&^«?=5)dS 

Jo 

= J  ~5f(dw®=5)d6 1  ‘(t=r+d)  o°(5rf=t)  />(S?  |  st°^  &  du°  =5)  dr  (34) 

The  first  equality  is  the  most  useful  one.  However  equation  (34)  is  useful  for  the  treatment  of 
TCGs  with  cycles  (see  subsection  4.4).  The  same  remark  applies  to  the  other  equations  involving 
P  distributions. 

2.2.  Cf:  lo  MET-BY  if 

P{C°&^\stQ  =  t&.duQ=d) 

=  f  "Sf  (if«f =5)  af  (^rf =c-d)P(S^\s/l=t-5  &du° = 5)  d5 

Jo 

= j  "5f  (dwf =5)  d5  J  ^p  '(t=r-5)  af(sff=T)  P{S°  I  ^rf =x  &  dii°  =5)  dz  (35) 

23.  Cf:  Iq  before  if 

/*(Cf  &  Sf  I  jfo  =  t  &  dwo = d) 

=J  5f(d«f=5)d6j  af(jrf=T)/’(5f|jrf=T&d«f=5)dT 


(36) 
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2.4.  C":  Iq  after  ij* 

/’(Cf&SfUrQ  =  t&duQ=d) 

=j^S^idu°=5)dS^'  ^a'^{st°=z)PiS^\st‘^=x&du°=d)dz  (37) 

2J.  Cf:  lo  DURING  if 

P(Cf&5f|«o  =  r&4/Uo=^ 

=  f  Sj'(duf=6)<i5r  a°(s/^=x)P{^\st°=x&du°=8)dx  (38) 

Jd  Jt+d-6 

2.6.  Cf:  lo  CONTAINS  if 

Pidl&^\sto  =  t&duQ=d)- 

=j^^(du°=5)dbj'  ‘‘  ^c!°(sr°=x)Pi:^\st°=x&.du°=5i)dx  (39) 

2.7.  Cf:  lo  STARTS  if 

P(d;&:^\sto  =  t&du^=d) 

=  f"5f(duf  =  5)o>f=0?(SfUff=f&duf=5)^i5 

•  (i 

=j^  5f(<i«f=5)<f6j  ^Kx=t)cs°(st1=x)P(^\st^=^&du°=5)dx  (40) 

2.8.  cf:  lo  STARTED-BY  if 

P{(fl&Sfl\st^  =  t&dUf,=d) 

=  f  ‘‘^{dul = S)  <7f(jrf = t)  F(5f  1 5rf =f &dMf = 5)  <^5 
Jo 

= J^‘'5f(^/uf  =5)^/5  j"’3‘(T=r)af(jrf=T)/»(Sf|.yff^&rfttf=5)^/T  (41) 

2.9.  cf:  lo  FINISHES  if 

P(Cf&Sfl5ro  =  r&c/Uo=^ 

=  f  “5f  (dwf = 5)  cf  (5  rf  =  r + d  -  5)  P(Sf  I  jff = f + d  -  5&  = 5)  d5 

•  if 

5f(duf=6)c^6j  p'('c=^+<3(-5)cf(5rf=T)F(Sf |jrf=T&^«f=6)dT  (42) 
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2.10.  C°;  FINISHED-BY  I® 

P{Cf^&^\stQ  =  t&duQ  =  d) 

=  f  5'^{du^  =  S)a°{st'^=t+d-5)P{^\st‘l=t+d-S&du^=8)db 
Jo 

J  ‘‘^(du^  =5)  j  "■p  ‘  (T=f-K/-5)  o°(st°=i)  P(5®  I  st°=^&du°=5)  dx  (43 ) 

2.11.  C®:  lo  OVERLAPS  ij* 

P((f,&.f,\st^  =  t&.du^=d) 

=  f**5®(du®=5)J5f'*^  <5%st^^^)P{fi\st°=x8cdu°=b)dx 

Jo  J  Max{t^*ii-6] 

=  f‘'5®(d«°=5)ii5f'^  a%t°=x)Pi^\s/l=x&dul=S)dx 
Jo  Jt*d-i 

■h 

^’’^{du°=S)dbj'*^a°{st°=x)P(f^\st°=x&du°=d^  (44) 

2.12.  Cf:  Iq  OVERLAPPED-BY  ij* 

P{C^&^\sto  =  t&dUg=d) 

= j"5®(i/M®=5)c/5j“"‘'’'*^“®  a®(jr°  =t)P(5®  |  sil=x&du° = 5)  dx 

=j^^idu°=8)dbj‘j5°ist°  =z)Pi^\st°=x&du^ = 6)i/t 
+ 

J  °°d°{dul=5)d5j  ^'^a®(jf°  =t)  P(5®  1  jr® = t  &  du^ = 5)  dx  (45) 

2.13.  C®:  Iq  equals  if 

P{C°Sc^\st^  =  t&du^=d) 

=  6%du!l=d)  o%st°=t)  P(S°  I  st°=t  &  du°=d) 

= j  "5®(^iu°=6)^i5  j  '~p2(5=4/,T=f)af(5?®=T)P(S®  1  ^r°=T  &  du^=S)dx  (46) 
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Appendix2:  Activity  Individual  Demand  Densities 


1.  Notations 

In  this  appendix  we  assume  that  a  posteriori  start  time  densities  have  already  been  computed  as 
described  in  section  4.  We  assume  that  none  of  these  densities  is  uniformly  zero,  otherwise  this 
would  indicate  unsatisfiability  of  the  current  CSP  and  the  incremental  scheduler  would  have  to 
backtrack.  As  already  mentioned  earlier,  a  posteriori  start  time  densities  can  then  be  normalized 
to  express  the  faa  that  each  activity  will  occur  once  (i.e.  each  activity  will  start  exactly  once). 
These  normalized  a  posteriori  densities  will  be  denoted* 

•  fixed-duration  activities: 

•  variable-duration  activities:  P f/^st=t&uiu~d&e.C ^ 

win  denote  the  a  priori  probability  that  uses  to  fulfill  its  resource  requirement 

Ry.  Djyy(t)  will  denote  A^’s  individual  demand  for  as  a  fimction  of  t.  This  is  the  probability 

that  A^  uses  at  time  t  to  fulfill  its  resource  requirement  R^.  The  computations  are  performed 
assuming  an  incremental  scheduler  whose  earlier  resource  reservations  are  non-preemptible. 
Therefore  the  demand  density  has  to  be  reshaped  so  that  it  does  not  overlap  with  earlier 
reservations.  We  propose  a  method  for  doing  so,  which  involves  two  steps. 

Finally  we  will  be  using  the  predicate  AVAIL(R4,y,t,t+dUjt)  which  returns  true  if  and  only  if 
resource  R^^  is  available  at  all  time  between  t  and  t+du^.  This  is  a  precondition  for  scheduling 
activity  A^  to  start  at  time  r,  if  A^^  is  to  use  resource  R^^. 

2.  Resource  Demand  Densities  Produced  by  Fixed-Duration  Activities 

The  probability  that  activity  Aj^  uses  Ri^^j  at  time  t  to  fulfill  its  resource  requirement  Rj^  is 
given  by  the  a  priori  probability  that  uses  Rj^j  to  fulfill  Rj^  multiplied  by  the  conditional 
probability  that  Aj.  is  active  at  time  r  given  that  it  uses  R^^j  to  fulfill  Ry.  It  turns  out  that  this 
latter  conditional  probability  may  be  uniformly  zero  for  some  resources  Rj^j  due  to  earlier 
reservations.  This  can  be  accounted  for  by  refining  the  a  priori  probabilities  (Rkij)-  D^,y(t)  is 
therefore  computed  in  two  steps: 

1.  In  the  first  step  we  compute: 

where  Pf/,sri^=x&Ci&...&CJRi^-=R^)  is  the  probability  that  st^=^  and  that  the 
temporal  relation  constraints  C,,...,C^  are  satisfied  given  the  activities’  a  priori  start 
time  distributions  and  given  that  R^^y  is  the  resource  used  to  fulfill  requirement  R^^^. 

This  probability  can  be  approximated  by  computing  I^’s  a  posteriori  start  time 
distribution  staning  from  an  a  priori  start  time  distribution  that  accounts  for  Ryy’s 
reservations.  We  do  so  by  replacing  (j^(sr^=r)  with  a^Csr^=r)xa[AVA/L(R^^j,/,r+du^)] 
in  the  computation  of  P(sq=t&C,&...&C„).  In  other  words: 

Pjv(srt-r&C,&...&C„IPi=P^-^=KP(sTt=x&C,&...&C„) 


a[AVAJL(R^,j,T,z-hdu^)] 


pa  5^ 
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where  K  is  a  normalization  factor^^. 

2.  In  the  second  step  the  a  priori  probabilities  Pr^  (R^y)  are  refined  to  account  for  the 
resource  reservations.  The  refined  probabilities  are  denoted  ).  Indeed,  due 

^ki 

to  earlier  reservations,  some  resources  Rj^y  have  a  posteriori  probabilities 
Pl^(st^=x&C^&...&CJR|^~R,^■p  that  are  uniformly  zero.  A^^’s  individual  demand  for 
these  resources  is  therefore  uniformly  zero  as  well.  Hence  one  can  use  the  new 
probabilities: 

p"‘P2(l?y  )z=  / 0  if  Pl^(st^=x&C^&....&CJR,=R^^J)  is  uniformly  0 
^  otherwise 

where  K^^  is  a  normalization  factor.  Notice  that,  for  each  R^,,  because  of  the 
consistency  checking  performed  after  the  computation  of  the  a  posteriori  start  time 
distributions,  we  are  guaranteed  at  this  point  to  have  at  least  one  resource  Rj,y  such 
that  P/^(st^=t<S:C Rfa=R*iy)  is  not  uniformly  zero.  One  can  then  compute: 

D~''w=pp^(Vj' 

practice  it  is  not  necessary  to  compute  one  can  just  compute  P^(st^=t&C,&...&C„,l 

,—^kip  check  if  it  is  uniformly  zero  or  not 

Finally,  notice  that  the  total  demand  is  given  by: 

(using  Fubini) 

=  (since  Pj,^  is  normalized) 

^Id 

and  hence  for  each  Ry  required  by  : 

£f“D;;^(TWT=<;u, 

I 

which  simply  expresses  that  an  activity  A^’s  total  demand  for  a  resource  Ry  is  equal  to  its 
duration  du;^.  This  duration  has  simply  been  distributed  over  time  and  over  several  resources 
(Ryy)  to  account  for  the  different  possible  schedules  of  the  aaivity. 
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Again  this  normalization  simply  expresses  that  the  activity  will  occur  exactly  once. 
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3.  Resource  Demand  Densities  Produced  by  Variable-Duration  Activities 

The  computations  in  the  case  of  variable-duration  activities  are  very  similar  to  the  ones  for 
fixed-duration  activities: 

1 .  In  the  first  step  one  computes  the  distributions 

a[AVA/L(/?,...T.T+e)] 

where  K  is  a  normalization  factor. 

2.  The  probabilities  Pn^CRfay)  are  refined  in  the  same  way  as  for  fixed-duration 
activities.  One  can  then  compute: 

Lastly,  using  Fubini’s  theorem,  one  can  check  that: 

Hence,  for  each  R^,-  required  by  A;^,  A^^’s  total  demand  is: 


which  is  A^’s  expected  duration  given  the  joint  start  time  and  duration  probability  density 

Acknowledgement 

We  like  to  thank  Nicola  Muscettola  and  Katia  Sycara  for  their  comments  on  this  work,  and  Joe 
Mattis  for  his  help  implementing  the  interface  to  the  scheduler. 


49 


References 


[Allen  83] 

J.F.Allen. 

Mainuining  Knowledge  about  Temporal  Intervals. 

Communications  of  the  ACM  26(1 1):832-843,  1983. 

[Allen  84] 

J.FAllen. 

Towards  a  General  Theory  of  Action  and  Time. 

Artificial  Intelligence  23(2):  123- 154, 1984. 

[Baker  74] 

K.R.  Baker. 

Introduction  to  Sequencing  and  Scheduling. 

WUey,  1974. 

[BeU  84] 

Colin  E.  Bell  and  Austin  Tate. 

Using  Temporal  Constraints  to  Restria  Search  in  a  Planner. 

Technical  Report  AIAI-TR-5,  Artificial  Intelligence  Applications  Institute, 
University  of  Edinburgh,  80  South  Bridge,  Edinburgh  EHl  IHN,  United 
Kingdom,  1984. 

[Dechter  88] 

Rina  Dechter  and  Judea  Pearl. 

Network-Based  Heuristics  for  Constraint  Satisfaction  Problem' 

Artificial  Intelligence  34(l):l-38, 1988. 

[Fox  83] 

M.  Fox. 

Constraint-Directed  Search:  A  Case  Study  of  Job-Shop  Scheduling. 

PhD  thesis,  Department  of  Computer  Science,  Camegie-Mellon  University, 
1983. 

[Freuder  82] 

E.C.  Freuder. 

A  Sufficient  Condition  for  Backtrack-fiee  Search. 

Journal  of  the  ACM  29(l):24-32,  1982. 

[Haralick  80] 

Robert  M.  Haralick  and  Gordon  L.  Elliott. 

Increasing  Tree  Search  Efficiency  for  Constraint  Satisfaction  Problems. 
Artificial  Intelligence  14(3):263-313,  1980. 

[Johnson  74] 

L.A.  Johnson  and  D.C.  Montgomery. 

Operations  Research  in  Production  Planning,  Scheduling,  and  Inventory 
Control. 

Wiley,  1974. 

[Lepage  78] 

G.  Peter  Lepage. 

A  New  Algorithm  for  Adaptive  Multidimensional  Integration. 

Journal  of  Computational  Physics  27:192-203,  1978. 

[LePape  87] 

Qaude  Le  Pape  and  Stephen  F.  Smith. 

Management  of  Temporal  Constraints  for  Factory  Scheduling. 

Technical  Report,  The  Robotics  Institute,  Carnegie  Mellon  University, 
Pittsburgh,  PA  15213, 1987. 

also  appeared  in  Proc.  Working  Conference  on  Temporal  Aspects  in 
Information  Systems,  Sponsored  by  AFCET  and  IFIP  Technical 
Committee  TC8,  North  Holland  Publishers,  Paris,  France,  May  1987. 

50 


[Mackworth  77] 

[Muscettola  87] 

[Nadel  86a] 

[Nadel  86b] 

[Nadel  86c] 

[Nudel  83] 

[Ow  84] 

[Purdom  83] 

[Rit  86] 

[Smith  83] 


A. K.  Mackworth. 

Consistency  in  Networks  of  Relations. 

Artificial  Intelligence  8(1):99-118,  1977. 

Nicola  Muscettola,  and  Stephen  Smith. 

A  Probabilistic  Framework  for  Resource-Constrained  Muti- Agent  Planning. 

In  Proceedings  of  the  Tenth  International  Conference  on  Artificial 
Intelligence,  pages  1063-1066.  1987. 

B. A.  Nadel. 

The  General  Consistent  Labeling  (or  Constraint  Satisfaction)  Problem. 

Technical  Report  DCS-TR-170,  Department  of  Computer  Science,  Laboratory 
for  Computer  Research,  Rutgers  University,  New  Brunswick,  NJ  08903, 
1986. 

B.A.  Nadel. 

Three  Constraint  Satisfaction  Algorithms  and  Their  Complexities:  Search- 
Order  Dependent  and  Effectively  Instance-specific  Results. 

Technical  Report  DCS-TR-171,  Department  of  Computer  Science,  Laboratory 
for  Computer  Research,  Rutgers  University,  New  Brunswick,  NJ  08903, 
1986. 

B.A.  Nadel. 

Theory-based  Search-order  Selection  for  Constraint  Satisfaction  Problems. 

Technical  Report  DCS-TR-183,  Department  of  Computer  Science,  Laboratory 
for  Computer  Research,  Rutgers  University,  New  Brunswick,  NJ  08903, 
1986. 

Bernard  Nudel. 

Consistent-Labeling  Problems  and  their  Algorithms:  Expected-Complexities 
and  Theory-Based  Heuristics. 

Artificial  Intelligence  21:135-178, 1983. 

Peng  Si  Ow. 

Heuristic  Knowledge  and  Search  for  Scheduling. 

Technical  Report,  Robotics  Institute,  Carnegie  Mellon  University,  Pittsburgh, 
PA  15213,  1984. 

(Ph.D.  Thesis). 

Paul  W.  Purdom,  Jr. 

Search  Rearrangement  Backtracking  and  Pol5momial  Average  Time. 

Artificial  Intelligence  21:117-133,  1983. 

Jean-Francois  Rit. 

Propagating  Temporal  Constraints  for  Scheduling. 

In  Proceedings  of  the  Sixth  National  Conference  on  Artificial  Intelligence, 
pages  383-388.  1986. 

Stephen  F.  Smith. 

Exploiting  Temporal  Knowledge  to  Organize  Constraints. 

Technical  Report,  Robotics  Institute,  Carnegie  Mellon  University,  Pittsburgh, 
PA  15213, 1983. 


51 


[Smith  85] 


[Smith  86] 


[Stone  86] 


[Stroud  71] 


[Thomas  83] 


[Tsang  87] 


Stephen  F.  Smith  and  Peng  Si  Ow. 

The  Use  of  Multiple  Problem  Decompositions  in  Time  Constrained  Planning 
Tasks. 

In  Proceedings  of  the  Ninth  International  Conference  on  Artificial 
Intelligence,  pzigts  1013-1015.  1985. 

S.  Smith,  M.  Fox,  and  P.S.  Ow. 

Constructing  and  Maintaining  Detailed  Produaion  Plans:  Investigations  into 
the  Development  of  Knowledge-Based  Factory  Scheduling  Systems. 

AI  Magazine  7(4):45-61,  Fall,  1986. 

Harold  S.  Stone  and  Paolo  Sipala. 

The  average  complexity  of  depth-first  search  with  backtracking  and  cutoff. 

IBM  Journal  of  Research  and  Development  30(3):242-258, 1986. 

A.  H.  Stroud. 

Series  in  Automatic  Computation:  Approximate  Calculations  of  Multiple 
Integrals. 

Prentice  Hall,  1971. 

Thomas  and  Finney. 

Calculus  and  Analytic  Geometry,  Sixth  Edition. 

Addison  Wesley,  1983. 

Edward  P.  K.  Tsang. 

The  Consistent  Labeling  Problem  in  Temporal  Reasoning. 

In  Proceedings  of  the  Seventh  National  Conference  on  Artificial  Intelligence. 
pages  251-255.  1987. 


[Valdes-Perez  87]  Raul  E.  Valdes-Perez. 

The  Satisfiability  of  Temporal  Constraint  Networks. 

In  Proceedings  of  the  Seventh  National  Conference  on  Artificial  Intelligence, 
pages  256-260.  1987. 

[Vere  83]  Stephen  Vere. 

Planning  in  Time:  Windows  and  Durations  for  Activities  and  Goals. 

IEEE  Transactions  on  Pattern  Analysis  and  Machine  Intelligence 
PAMI-5"(3):246-267, 1983. 

[Vilain  86]  Marc  Vilain  and  Heruy  Kautz. 

Constraint  Propagation  Algorithms  for  Temporal  Reasoning. 

In  Proceedings  of  the  Sixth  National  Conference  on  Artificial  Intelligence, 
pages  377-382.  1986. 


